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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2635 


AN  ANALYSIS  OF  LAMINAR  FREE  -CONVECTION  FLOW  AND  HEAT  TRANSFER 
ABOUT  A  FLAT  PLATE  PARALLEL  TO  THE  DIRECTION 
OF  THE  GENERATING  BODY  FORCE 


By  Simon  Ostrach 


SUMMARY 

Tlie  free -convection  flow  and  heat  transfer  (generated  by  a  "body 
force)  about  a  flat  plate  parallel  to  the  direction  of  the  body  force 
are  formally  analyzed  and  the  type  of  flow  is  found  to  be. dependent  on 
the  Grashof  number  alone.  For  large  Grashof  numbers  (which  are  of 
interest  in  aeronautics),  the  flow  is  of  the  boundary -layer  type  and  the 
problem  is  reduced  in  a  formal  manner,  which  is  analogous  to  Prandtl' s 
^  forced -flow  boundary -layer  theory,  to  the  simultaneous  solution  of  two 

ordinary  differential  equations  subject  to  the  proper  boundary  conditions. 

S  Velocity  and  temperature  distributions  for  Prandtl  numbers  of  0.01, 

0.72,  0.733,  1,  2,  10,  100,  and  1000  are  computed  and  it  is  shown  that 
velocities  and  Nusselt  numbers  of  the  order  of  magnitude  of  those 
encountered  in  forced -convection  flows  may  be  obtained  in  free -convection 
flows.  The  theoretical  and  experimental  velocity  and  temperature  distri¬ 
butions  are  in  good  agreement. 

A  flow  and  a  heat-transfer  parameter,  from  which  the  important  phys¬ 
ical  quantities  such  as  shear  stress  and  heat-transfer  rate  can  be  com¬ 
puted,  are  derived  as  functions  of  Prandtl  number  alone.  Comparison  of 
theoretically  computed  values  of  the  heat -transfer  parameter  with  values 
obtained  from  an  approximate  calculation  and  experiments  yielded  good 
agreement  over  a  large  range  of  Prandtl  number.  Agreement  between  the 
theoretical  values  and  those  obtained  from  a  frequently  used  semiempiri- 
cal  heat-transfer  law  was  good  only  in  restricted  Prandtl  number  ranges 
(depending  on  an  arbitrary  constant). 


INTRODUCTION 

Two  important  types  of  fluid  flow  problems  involving  heat  transfer 
are  those  of  forced  and  those  of  free  convection.  By  forced -convection 
flow  is  meant  flows  maintained  mechanically  as,  for  example,  by  a  pres¬ 
sure  drop  or  an  agitator.  Free -convection  flow,  on  the  other  hand, 
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results  from  the  action  of  "body  forces  on  the  fluid,  that  is,  forces 
which  are  proportional  to  the  mass  or  the  density  of  the  fluid.  The 
flow  is  generally  produced  in  the  following  manner;  Consider,  for 
example,  a  fixed  object  (such  as  a  plate)  in  a  quiescent  fluid  subject 
to  a  body  force.  When  the  plate  is  at  the  same  temperature  as  the  sur¬ 
rounding  fluid,  the  body  forces  acting  on  the  fluid  are  in  equilibrium 
with  the  hydrostatic  pressure  and  no  flow  ensues  in  the  steady  state. 

If  a  temperature  gradient  normal  to  the  body  force  is  imposed  by  heating 
(or  cooling)  the  plate,  there  will  exist  a  defect  (or  excess)  of  body 
force  because  of  the  decreased  (or  increased)  density,  with  the  fluid 
closer  to  the  plate  having  the  greater  defect  (or  excess)  than  that  away 
from  the  plate.  This  unbalance  of  the  forces  causes  the  fluid  to  be 
accelerated  with  the  particles  nearer  the  plate  moving  more  rapidly  than 
those  farther  from  the  plate.  Free -convection  flow  has  usually  been  con¬ 
sidered  to  be  generated  in  a  gravitational  field  where  the  previously 
mentioned  defect  or  excess  of  body  force  was  the  Archemedian  (buoyancy) 
force.  However,  since  centrifugal  forces  are  also  proportional  to  the 
fluid  density,  free -convection  flows  can  also  be  set  up  by  the  action  of 
such  forces.  (See  reference  1  for  a  more  explicit  discussion  of  the 
development  of  free-convection  flows  by  centrifugal  forces.) 

Free -convection  flows  produced  by  centrifugal  forces  are  now  of 
practical  importance  in  aeronautics  because  many  aircraft  propulsion  g 

systems  contain  components  (such  as  gas  turbines  and  helicopter  ram  jets) 
which  rotate  at  high  speeds  and  in  which  heat  is  being  transferred.  The 
method  of  free-convection  cooling  of  gas-turbine  rotor  blades  where  the 
centrifugal  forces  create  a  free-convection  flow  of  the  coolant  in  the 
blade  passages  is  an  example  of  a  practical  application  of  the  free- 
convection  phenomenon  in  aeronautics.  Also,  free-convection  flow  due 
to  centrifugal  force  is  superimposed  on  the  flow  through  helicopter  ram 
jets  and  on  the  flow  of  cooling  air  in  hollow  rotor  blades  of  air-cooled 
turbines  and,  under  proper  conditions,  can  appreciably  influence  the 
resultant  flow  and  heat  transfer. 

As  a  simplification  of  the  many  free-convection  problems  which  are 
now  of  some  consequence  in  aeronautics,  consideration  is  here  given  to 
the  special  case  of  free-convection  flow  about  a  flat  plate  parallel  to 
the  direction  of  the  generating  body  force.  The  experimental  and  theo¬ 
retical  considerations  of  Schmidt  and  Beckmann  (reference  2)  concerning 
the  free-convection  flow  of  air  subject  to  the  gravitational  force  about 
a  vertical  flat  plate  constitute  the  most  complete  treatment  of  this 
subject  up  to  the  present  time.  Eckert  (reference  l)  as  well  as  others 
has  further  verified  and  extended  the  experimental  results  of  Schmidt 
and  Beckmann,  and  Schuh  (reference  3)  has  extended  the  numerical  calcu¬ 
lations  by  computing  the  velocity  and  temperature  distributions  fbr 
several  Prandtl  numbers  different  from  that  for  air.  However,  all  the 
theoretical  work  in  these  references  is  based  on  the  incompressible 
equations  in  which  the  density  (or  temperature)  variation  is  introduced  1 
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f  i-n  the  "buoyancy  term  alone.  Yarious  terms  are  omitted  from  the  equations 

at  the  start  on  the  basis  of  either  intuitive  arguments  or  no  arguments 
at  all.  Although  a  theoretical  development  made  in  such  a  manner  led  to 
good  final  results,  the  significance  of  all  the  important  factors  assoc¬ 
iated  with  the  free -convection  flow  phenomenon  is  not  obtained  from  such 
an  analysis. 

The  problem  of  free-convection  flow  as  produced  by  a  body  force 
about  a  flat  plate  in  the  direction  of  the  body  force  was  studied  at 
the  NACA  Lewis  laboratory  and  is  treated  in  a  formal  and  more  general 
manner  herein.  The  method  used  is  somewhat  similar  to  that  used  in 
reference  4  wherein  Consideration  was  given  to  the  free-convection  flow 
at  high  Grashof  numbers  in  a  horizontal  cylinder  which  had  a  variable 
surface-temperature  distribution.  The  application  of  this  method  to 
the  present  problem  leads  to  a  development  which  is  analogous  to  Prandtl's 
treatment  of  high  Reynolds  number  forced -convection  flows.  Although  the 
final  equations  obtained  by  this  method  are  the  same  as  those  of  Schmidt 
and  Beckmann,  this  more  general  approach  not  only  clearly  demonstrates 
the  significance  of  all  the  important  parameters  and  assumptions  and 
hence  leads  to  a  better  understanding  of  this  type  of  flow  but  also 
indicates  the  quantitative  limitations  of  the  theory.  In  addition,  the 
numerical  solutions  of  references  2  and  3  are  herein  extended  to  cover 
*  a  more  complete  range  of  variables .  The  new  calculations  yield  infor¬ 

mation  on  the  free-convection  flow  for  Prandtl  numbers  corresponding  to 
those  of  liquid  metals,  gases,  liquids,  and  very  viscous  fluids. 


ANALYSIS 

Statement  of  Problem  and  Basic  Equations 

The  steady -state  equations  expressing  the  conservation  of  mass, 
momentum,  and  energy  for  a  compressible,  viscous,  and  heat-conducting 
fluid  subject  to  a  body  force  together  with  an  equation-  of  state  govern 
the  flow  and  associated  temperature  distribution  about  the  plate.  These 
equations  in  Cartesian  tensor  notation  are  (see  reference  5),  respec¬ 
tively, 

jL  (pu.)  =  0  (1) 

Sx .  J 


* 
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M  2 

\Slj  +  dxlyl  "  3  W  _ 

(3) 

P  =  P(P,T)  (4) 

(A  complete  list  of  the  symbols  used,  herein  is  given  in  appendix  A.) 

For  the  two -dimensional  case,  equations  (l)  to  (4)  represent  a  system 
of  five  equations  in  the  five  dependent  variables  U-j_;  Ug,  P,  P,  and  T. 
For  later  use,  equation  (4)  can  he  written 

dp  =  p(K  dP  -  (3  dT)  (4a) 

where  K  and  (3  are  the  coefficients  of  isothermal  compressibility  and 
volumetric  expansion,  respectively  (see  reference  6) .  In  addition  to  a 
general  state  equation,  such  as  is  given  in  equations  (4)  or  (4a),  it 
will  be  convenient  at  times  in  the  discussion  to  refer  to  some  specific 
state  equation.  To  this  end,  the  equation  of  state  for  an  ideal  gas 

P  =  pRT  (4b) 


pcpuj 


Sx 


=  u. 


dX 


ax. 


dT 

dx 


+  M- 


will  be  used. 

Particular  consideration  is  here  given  to  the  two-dimensional  free- 
convection  flow  about  a  semi-infinite  vertical  flat  plate.  The  X-^axis 
of  the  coordinate  system  is  taken  along  the  plate  and  the  Xg-axis,  nor¬ 
mal  to  it.  No  distinction  is  made  as  to  the  specific  type  of  body  force 
acting,  for  example,  gravitational  or  centrifugal,  but  the  force  is 
assumed  to  be  acting  in  the  vertical  direction  only  (that  is,  parallel 
to  the  plate) .  Centrifugal  and  Coriolis  forces  which  are  connected  with 
flows  on  curved  paths  and  with  rotating  systems  generally  vary  with 
position  and  velocity.  However,  in  order  not  to  make  the  analysis 
unduly  complicated,  the  body  force  is-  taken  to  be  constant. 

In  order  to  define  the  problem  clearly,  a  choice  must  still  be  made 
of  the  position  of  the  origin  of  the  coordinate  system.  Before  making  a 
definite  decision  on  this  point,  note  that  for  constant  plate  tempera¬ 
tures  there  are  four  permutations  of  the  body -force  direction  (either 
upward  or  downward)  and  the  plate  thermal  condition  (either  heated  or 
cooled)  which  will  lead  to  free -convection  flows.  Once  the  position  of 
.the  edge  of  the  plate,  which  is  also  to  be  the  origin  of  the  coordinate 
system,  is  decided,  there  are  two  combinations  of  the  body -force  direc¬ 
tion  and  plate  thermal  condition  that  will  yield  flows  which  proceed 
away  from  the  edge .  It  is  this  type  of  flow  that  is  amenable  to  the  type 
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of  analysis  to  be  made  here.  This  point  will  he  more  fully  discussed 
subsequently .  If  the  edge  of  the  plate  (recall  that  a  semi -inf inite 
plate  has  but  one  edge)  is  taken  at  the  bottom  of  the  plate  (that  is, 
the  plate  extends  to  +®  in  the  X-^ -direction) ,  the  two  combinations  lead¬ 
ing  to  flows  in  the  proper  direction  (upward  in  this  case)  are,  respec¬ 
tively,  the  body  force  acting  downward  with  a  heated  plate  and  the  body 
force  acting  upward  with  a  cooled  plate.  The  equations  developed  for  one 
of  the  cases  reduce  directly  to  those  for  the  other.  The  remaining  two  - 
permutations,  namely,  the  body  force  acting  downward  with  a  cooled  plate 
and  the  body  force  acting  upward  with  a  heated  plate,  would  yield  flows 
which  proceed  downward  or  toward  the  edge  of  the  plate  if  this  edge  were 
taken  at  the  bottom  pf  the  plate.  This  type  of  flow  would  violate  a 
physical  condition  of  the  problem  which  states  that  the  flow  starts  at 
the  plate  edge.  The  latter  combinations  hence  will  not  be  considered 
further. 

Because  the  two  acceptable  configurations  can  be  reduced  essentially 
to  one,  for  the  development  to  be  given  here,  the  origin  of  the  coordinate 
system  will  be  taken  at  the  bottom  of  a  heated  plate,  with  the  body  force 
acting  downward.  The  assumption  is  now  made  that  the  viscosity  and  heat- 
conductivity  coefficients  are  functions  of  the  temperature  alone  and  obey 
the  following  laws: 

>-(£)"  -40T  <s> 

The  choice  of  the  body -force  direction  together  with  equations  (5)  alters 
equations  (2)  and  (3)  so  that  they  become 


(7) 

Note  that  the  only  nonzero  component  of  the  body  force  is  the  X^-component 


Boundary  Conditions 

The  boundary  conditions  associated  with  the  given  problem  are  that: 
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(a)  The  fluid,  must  adhere  to  the  plate  (the  no-slip  condition  of  tr 

viscous  flows)  and  the  plate  must  he  a  streamline,  or  mathematically, 

UiO^O)  =  U2(I1,0)  =  0  (8) 

(h)  The  temperature  of  the  fluid  at  the  plate  must  he  equal  to  the 
plate  temperature ;  that  is, 

T(X1;0)  =  T0  (9) 

(c)  The  velocity  U]_  at  large  distances  from  the  plate  must  he 
undisturbed ,  or 

Uxdi,-)  =  0  (10) 

(d)  The  temperature  at  large  distances  from  the  plate  must  he  equal 
to  the  undisturbed  fluid  temperature }  or 

T  (*!,-)  =  T^ 

Simplification  of  Equations 

Let  a  small  quantity  e  no w  he  defined  as 

e  =  P(T0-T.)  (12) 

which  is  a  measure  of  the  magnitude  of  temperature  variation  in  the  flow 
field.  The  coefficient  of  volumetric  expansion  (3  is  generally  of  the 
order  of  magnitude  between  10”^  and  10”^  (see  table  15  of  reference  7 , 
for  example)  and  for  gases,  p  =  l/T.  (Thus,  for  gases,  if  p  is  taken 
to  he  constant,  e  =  (Tq-T«>)/T»> ;  that  is,  e  is  the  relative  temperature 

difference.)  The  coefficient  p  will  he  assumed  constant.  Because  in 
the  steady  state  flow  ensues  only  when  there  is  a  temperature  variation 
in  the  fluid,  the  free -convection  velocity  should  then  depend  directly 
on  e,  and  the  variations  in  pressure  and  density  (from  the  static, 

6=0,  case)  due  to  the  temperature  differences  should  also  depend  on  6 
Thus 

Ui  =  6 

P  =  Ps 
P  =  Ps 

T  =  To, 


(13) 

+  T®6CT 

(14) 

+  P«6<P 

(15) 

( 1  +  cd) 

(16) 

(11) 
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where  -fx  denotes  the  X1-component  of  the  tody  force  per  unit  mass, 
ui,a,  q>,  and  B  denote  dimensionless  functions  (which,  in  general,  can 
he  functions  of  e),  2  is  some  characteristic  length  (for  example,  the 
distance  from  the  edge  of  the  plate  to  the  point  of  interest),  PQ  and 
Ps  jare  the  pressure  and  the  density,  respectively,  for  the  static  case, 
(Uj_  =  0  or  e  =  0),  and  Pw  and  denote  constant  values  of  the 

pressure  and  the  density  (that  is,  the  values  if  no  force  field  were  pre¬ 
sent)  defined  by  the  state  equation  (in  the  case  of  a  gas,  in  particular 
P®  =  ftoEToo) .  Because  there  is  no  characteristic  velocity  associated  with 

the  type  of  flow  under  consideration,  the  velocity  is  dimens ionalized  by 
the  factor  given  in  parentheses  on  the  right  side  of  equation  (13).  This 
specific  factor  is  chosen  because  it  leads  to  simpler  equations. 

In  order  to  determine  the  static  quantities,  it  will  at  first  be 
convenient  to  consider  the  particular  case  of  a  gas.  The  problem  is 
then  considered  with  the  temperature  uniform  throughout  the  flow  field 
at  the  value  T„  (therefore  there  will  be  no  flow  and  U.  =  0).  For 
this  situation,  equations  (4b)  and  (6)  become 

Pg  ss  Pg  RT«o  (17) 

and 


Sx- 


+  ps  fX  =  0 


bx. 


=  0 


> 


J 


(18) 


(it  should  be  noted  that  equation  (18)  expresses  the  physical  fact  pre¬ 
viously  stated  that  the  body  force  and  hydrostatic  pressure  are  in  equil¬ 
ibrium  for  the  static  case.)  Substitution  of  equation  (17)  into  equa¬ 
tion  (18)  leads  to 


Ps  =  So  exp 


(19) 


and  equation  (19)  together  with  equation  (17)  and  the  equation  defining 
Poo  and  po,  yields 


Pg  =  fto  exp 


Poo  ©xp 


(20) 


If  the.  exponential  in  equation  (20)  is  expressed  in  terms  of  its  series 
expansion,  that  equation  becomes 
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(21) 


A  computation  of  the  second  term  in  the  brackets  of  equation  (21)  for  the 
case  of  air  under  normal  conditions  with  f  x  =  g  and  the  fact  that  Xj_ 
is  of  the  order  of  magnitude  Z  show  that  p^gZ /~Pn  ~  10”^z/f oot .  For 
the  type  of  problem  under  consideration,  Z  will  always  be  of  unit  order 
of  magnitude  so  that  even  if  the  body  force  fx  represents  a  centrifugal 
force  many  times  that  of  gravity,  the  inequality  p*  f  jXq/Pa>  «  1  may 
still  be  satisfied.  Thus,  in  the  subsequent  development  it  will  be 
assumed  that  ps  s  p,,,.  This  assumption  which  was  justified  by  the  com¬ 
putation  for  the  case  of  a  gas  is  expected  to  be  reasonable  for  other 
fluids  as  well.  The  physical  interpretation  of  this  assumption  is  that 
under  static  conditions  (e=  0),  the  density  (or  pressure)  is  not 
affected  by  the  force  field. 


In  order  that  all  quantities  in  the  following  equations  be  dimen¬ 
sionless,  it  is  further  necessary  to  define  x.^  =  X-^/Z,  where  the  x^^ 
axe  now  dimensionless  space  coordinates.  Substituting  these  new  coor¬ 
dinates  along  with  equations  (13)  to  (16)  into  equations  (l),  (6),  (7), 
and  (4a)  and  noting  equations  (18)  and  that  pg =  p„  yield,  on  neglec- 
tion  of  terms  of  higher  order  in  e  compared  with  those  of  order  c, 

du]_ 

SX]_ 


dug 

dx2 


=  0 


(22) 


/  -  du;L  du-A 

8r  57  *  -■  4ui  ■ 

ui  +  uz  sy  ■  4u2 


-  KGr 


axq 

da 


ax. 


(23) 

(24) 


Gr  *  (ui  ^ + u2  19 = ^ Gr  *  (ui  + us  ly + Ae  (25) 


d<p  =  K  I^da  -  PT*.  d© 


(26) 


where  Poo/ftofyZ  =  NGr  and  the  Grashof  number  Gr  and  the  Prandtl  number 
Pr  are  defined  as 


p„2fxZ3e 


Pr  = 


C  LL 
P'  co 


V 


l 


» 


r 


Gr  = 


and 


r 
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Physically  speaking,  the  Grashof  number  represents  the  ratio  of  the  body 
forces  to  the  viscous  forces.  The  free-convection  flows  of  interest  here 
are  those  associated  with  large  Grashof  numbers .  The  factors  K  and  |3 
in  equation  (26)  may  well  be  taken  to  be  constants  (see  reference  6). 

The  boundary  conditions  (equations  (8)  to  (ll))  in  terms  of  the  new 
dimensionless  variables  are 


u1(x1, 0)  =  U2(X]_,0)  =  0 

(27) 

e(li>0)  *  sfc 

(28) 

ul(xl>“)  =  o 

(29) 

e(*L,-)  =  0  (30) 

Thus,  to  a  first  approximation,  equations  (22)  to  (26)  together  with  the 
houndary  conditions  replace  the  original  equations  and  houndary  conditions 
(Note  that  for  gases,  PT^  =  1.) 

The  prime  assumption  that  has  been  made  in  this  analysis  is  that  the 
higher-order  terms  in  e  were  negligible,  which  implies  that  e  is 
small  and,  consequently,  that  the  temperature  difference  or  p  is 
moderately  small.  It  is  a  consequence  of  this  assumption  alone  that  the 
basic  equations  were  simplified  to  equations  (22)  to  (26)  wherein  the 
viscosity  term  in  the  energy  equation  is  neglected  and  the  only  coupling 
of  the  momentum  and  energy  equations  occurs  by  means  of  the  body-force 
term  in  equation  (23).  As  a  result  of  this  assumption,  the  variations 
of  the  viscosity  and  heat -conductivity  coefficients  with  temperature  are 
also  negligible.  Without  any  discussion,  the  authors  of  reference  2 
start  directly  from  simplified  equations  of  the  same  form  wherein  the 
pressure  terms  In  the  energy  equation  were  also  neglected.  In  refer¬ 
ence  3  some  intuitive  arguments  are  given  to  justify  the  simplified 
equations . 


It  is  now  convenient  to  revert  to  the  more  familiar  notation  where 
x  =  x-j_,  y  =  x2,  u  as  U]_ ,  and  v  =  u2.  Equation  (22)  implies  the  existence 
of  a  stream  function  ¥  such  that 


u 


and 


v 


-*x  (3D 


where  subscripts  denote  differentiation.  Applying  equation  (31)  to 
equations  (23)  to  (25),  yields,  respectively, 
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Gr  ^  *yYxy  '  *xYyy  +  N°x  (32^ 

Gr  (A*x)  ~  “^y^xx  +  ^x^xy  +  ^y  (33) 

(34) 


Yy(x,  0)  =  Yx(x,0)  =  0 

(35) 

0(x,O)  =  — i— 

(36) 

Yy(x,«)  =  e(x,»)  =  o 

(37) 

Equations  (32)  to  (34)  and.  equation  (26)  form  the  system  of  equations 
for  the  four  unknown  functions  Y,  0,  <p,  and  a  of  the  problem.  The 
system  is  nonlinear,  and  therefore  further  simplification  of  the  equations 
would  he  desirable.  Just  as  in  the  case  of  forced -convection  flows  where 
the  Reynolds >  number  determines  the  type  of  flow  or,  in  mathematical  terms, 
the  type  of  solution,  the  Grashof  number  is  the  prime  factor  for  free- 
convection  flows.  For  the  case  of  small  Grashof  number,  it  can  be  seen 
from  equations  (32)  to  (34)  that  a  perturbation  in  the  small  parameter 
Gr  will  yield  a  system  of  linear  equations .  For  Grashof  numbers  of  unit 
order  of  magnitude,  no  further  important  simplification  can  be  made  and 
the  solutions  would  have  to  be  obtained  numerically.  For  the  other  lim¬ 
iting  case,  that  of  large  Grashof  numbers  (which  is  the  case  under  con¬ 
sideration  herein),  it  would,  at  first  thought,  appear  that  some  simpli¬ 
fication  could  be  obtained  by  performing  a  perturbation  in  the  small 
parameter  l/Gr.  However,  this  would  then  imply  that  the  term  containing 
the  highest-order  derivatives  (the  left  term  in  equations  (32)  to  (34)) 
could,  among  others,  be  neglected.  (This  argument  would  also  imply  that 
the  body -force  term  <p  in  equation  (32),  which  is  essentially  causing 
the  flow,  could  also  be  neglected.)  The  omission  of  the  highest-order 
derivatives  from  consideration,  however,  would  lead  to  solutions  which 
would  not  satisfy  all  the  boundary  conditions.  Problems  of  this  type 
are  referred  to  as  singular  perturbation  problems.  For  further  dis¬ 
cussions  of  singular  perturbation  problems  see  references  8  and  9. 

Equations  in  which  a  small  parameter  multiplies  the  highest-order 
terms  are  said  to  be  of  the  boundary -layer  type,  because  in  order  for 
solutions  which  satisfy  all  the  boundary  conditions  to  be  obtained, 
the  highest-order  terms  must  be  considered  near  the  boundary.  This 
fact  implies  the  existence  of  a  thin  region,  called  the  boundary  layer, 
wherein  the  functions  vary  rapidly  from  the  value  at  the  boundary  to 
that  in  the  flow  outside  this  layer.  The  conclusion  to  be  drawn  from 


WW  =  V*  -  *x0y  -  yfr  (¥ yff*  -  W 

The  boundary  conditions  (equations  (27)  to  (30))  become 
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the  preceding  discussion  is  that  for  large  Grashof  numbers  the  f low 
is  of  the  boundary-layer  type.  Schmidt  and  Beckmann  (reference  2) 
also  made  the  boundary-layer  assumptions  in  their  theoretical  develop¬ 
ment,  and  these  assumptions  were  justified  on  the  basis  of  their  experi¬ 
mental  observations.  The  Grashof  numbers  for  their  experiments  were  of 
the  order  of  8xl06 . 

In  view  of  the  fact,  previously  discussed,  that  highest-order  deri¬ 
vatives  of  each  dependent  variable  as  well  as  of  those  terms  of  physical 
importance  (as,  for  example,  the  body -force  term)  must  be  retained  in  the 
boundary  layer,  it  is  convenient  to  make  both  sides  of  each  of  the  equa¬ 
tions  of  the  same  order  in  Gr.  In  this  way,  as  will  be  shown,  the 
equations  will  be  simplified  further.  It  is  thus  convenient  to  make  the 
following  transformations  in  the  system  of  equations  (32)  to  (34)  and 
(26)  and  then  to  retain  only  the  dominant  parts  (that  is,  those  multi¬ 
plied  by  Gr  to  the  highest  power)  of  each  individual  term. 

Let  y  =  Gr*1  y,  Y  =  Gr3  f ,  a  =  Gr*3  a,  <p  =  <p ,  and  0  =  0.  Then 
equations  (32)  to  (34)  and  (26)  become 


(38) 


(39) 


orS+3r~lfOT  - <Jr'1 5  -  Gr2B+2r<¥*?  -  +  H  0rts 


-Ors+2r-1  T  xyy  =  Or2B+r  (-t#„  +  +  N  Or**1?,-; 


J^xx  1  xy- 


% 


G£ 

Pr 


2r-l 


GrS+r  -  V,) 


7  Grt+a+r 
7-1 


(T  y5X  -  T2CTy)  (40) 


d<p  =  EP^Gr^  dcr  -  pT„.d0  (41) 

It  can  now  be  seen  that  by  proper  choice  of  r,  s,  and  t  a  transforma¬ 
tion  of  the  type  given  provides  a  means  for  making  the  important  terms 
in  the  differential  equations  of  the  same  order  in  Gr.  Thus  if  r  =  l/4, 
s  =  -3/4,  and  t  =  -1,  equations  (38)  to  (41)  become 


Y-Y  - 

IyI  xy - 


YxYyy 


+  No. 


(42) 


No_ 

y 


_i 

=  0(Gr  2)  £  0 


(43) 


Tr  Vi 


dq>  +  pT^d©  =  0 


(44) 


(45) 
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More  generally,  if  N  is  very  much  different  frcm.  unit  order  of  magni¬ 
tude,  a  value  of  t  can  always  he  chosen  (depending  on  N)  such  that 
equations  (42)  to  (45)  are  obtained.  (For  any  negative  t  less  than  -1, 
the  last  term  of  equation  (42)  -will  also  disappear.) 

There  are  now  several  important  points  to  he  .discussed  concerning 
the  transformation  just  made  and  the  resulting  simplified  equations. 
First,  it  should  he  noted  that  the  transformation  is  merely  a  formal 
expression  of  the  boundary -layer  assumptions  first  made  by  Prandtl  and 
hence  the  solutions  will  he  asymptotic  for  large  Gr.  Second,  the  second 
equation  of  motion  here  also  reduces  to  state  that  the  pressure  across 
the  boundary  layer  is  constant.  Third,  the  pressure  terms  in  the  energy 


and  state  equations  are  here  found  to  be  negligible, 
a  priori  assumptions  made  by  others  from  the  physics 
Finally,  note  that  integration  of  the  general  state 
of  pressure)  as  now  given  by  equation  (45)  leads  to 

This  fact  verifies 
of  the  problem, 
equation  (independent 

cp  +  pT^e  =  0 

(46) 

where  the  constant  of  integration  has  been  taken  as 
of  generality.  For  the  particular  case  of  a  gas,  p 
tion  (46)  becomes 

zero  without  any  loss 
=  l/Tm  so  that  equa- 

q>  +  0  =  0 

The  boundary  conditions  (equations  (35)  to  (37))  now  can  be  written 

fy(x,0)  =  Tx(x,0)  =  0 

(47) 

(48) 

T_(x,-)  =  =  0 

(49) 

If  now  it  is  assumed  that  <jx  =  0  in  equation  (42)  since  consideration 

is  here  being  given  to  a  flat  plate,  and  if  <p  is  eliminated  from  equa¬ 
tion  (42)  "by  use  of  equation  (46),  there  results  the  system  of  equations 

%y  +  PT.e  =  If  Xy  -  *x*yy 

(50) 

0-y  =  Pr(Ty0x  -  Tx0y) 

(51) 

Thus  the  problem  has  been  reduced  to  the  solution  of  the  two  simultaneous 
partial  differential  equations  (equations  (50)  and  (51))  subject  to  the 
boundary  conditions  (equations  (47)  to  (49)). 
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Final  simplification  of  the  equations  is  made  by  application  of  the 
so-called  similarity  transformation  of  boundary- layer  theory.  Thus,  let 

n  =  —z-i  (52) 

(4z)^ 

and 

3 

?  =  (4x)4  P(ti)  (53) 


9  -  '  <»> 

Then  equations  (50)  and  (51)  are  reduced  to  the  following  ordinary  dif¬ 
ferential  equations: 

F"  ’  +  3FF’  *  -  2F’2  +  H  s  0  (55) 

H"  +  2Pr  F  H’  =0  (56) 

where  the  primes  denote  differentiations  with  respect  to  i\.  The  reci¬ 
procal  one -fourth  power  similarity  as  given  in  equation  (52)  is  charac¬ 
teristic  of  free -convection  flows  just  as  the  reciprocal  square-root 
type’  is  characteristic  of  the  forced-convection  flows.  The  boundary 
conditions  became 

F’ (0)  =  F(0)  =  0  (57) 

H(0)  =1  (58) 

F’  (•)  =  H(«)  =  0  (59) 

The  use  of  a  transformation  like  equation  (52)  essentially  specifies  an 
additional  boundary  condition,  namely,  that  the  conditions  to  be  satis¬ 
fied  at  y  =  »(orTj!=«»)  should  also  be  satisfied  at  x  =  0.  It  is 
for  this  reason  that  the  flows  previously  discussed  which  would  flow 
toward  the  edge  (downward)  are  not  amenable  to  this  type  of  analysis, 
for  such  flows  would  violate  this  additional  condition,  which  essentially 
states  that  the  boundary -layer  development  starts  at  the  edge  of  the 
plate. 


Solution  of  the  Boundary -Value  Eroblem 

The  solutions  of  the  simplified  equations  (55)  and  (56),  satisfying 
the  boundary  conditions  as  given  by  equations  (57)  to  (59),  were  obtained 


14 


MCA  TN  2635 


ly  use  of  an  IBM  Card  Programmed  Electronic  Calculator.  A  detailed 
account  of  the  procedure  followed  in  the  determination  of  the  unknown 
functions  is  presented  in  appendix  B  by  Dr.  Lynn  U.  Albers.  The 
functions  F  and  H  together  with  their  derivatives  are  given  in 
table  I  for  Prandtl  numbers  of  0.01,  0.72,  0.733,  1,  2,  10,  100,  and 
1000.  Even  though  the  Prandtl  number  for  air  is  taken  as  0.72  in  this 
report,  the  solutions  for  Pr  =  0.733  were  also  computed  and  are  presented 
as  a  check  with  the  Schmidt -Beckmann  calculations  wherein  the  value  of 
Prandtl  number  of  0.733  was  used.  The  particular  values  of  the  Prandtl 
numbers  given  were  chosen  to  correspond  to  those  for  liquid  metals,  gases, 
liquids  (such  as  water  and  oil),  and  very  viscous  liquids  (such  as  glycerin 
or  oils  at  very  low  temperatures ) . 


RESULTS 

Velocity  and  Temperature  Distributions 

By  means  of  the  various  transformations  made  in  the  analysis  it  can 
easily  be  verified  that 


and 


where 


U 


TJX 

v_ 


2Vp(T0-Tlo)fzX  2-VGrJ 


=  F'(r}) 


n 


_.  i 

fx(Tp -Tg>)P  4 

4v„2X 


(60) 


(61) 


(62) 


Equations  (60)  to  (62)  relate  the  physical  quantities  to  the  dimension¬ 
less  functions  F  and  H  which  are  now  known.  The  dimensionless  velo¬ 
city  and  temperature  distributions  as  given  by  equations  (60)  and  (61) 
are  presented  in  figures  1  and  2,  respectively,  as  functions  of  rj  for 
the  various  values  of  Prandtl  number.  The  computations  made  here  agree 
with  those  for  Pr  =  0.733  as  given  in  reference  2  up  to  the  third  sig¬ 
nificant  figure.  For  Pr  =  10,  100,  and  1000,  the  present  results  agree 
in  general  with  those  of  reference  3 .  Since  only  curves  are  presented 
in  reference  3,  the  precision  of  the  agreement  cannot  be  stated. 


t 


t 
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The  maximum  values  of  the  dimensionless  velocity  distributions  occur 
at  larger  values  of  the  argument  tj  as  the  Prandtl  number  decreases  and 
the  velocities  decrease  with  increasing  Pr.  It  should  also  he  noted 
that  the  dynamic  and  thermal  boundary -layer  thicknesses  can  be  estimated 
from  the  abscissas  of  figures  1  and  2,  respectively,  and  that  for  Pr»  1 
the  velocity  boundary  layer  is  much  thicker  than  the  thermal  boundary 
layer . 

The  occurrence  of  fx  (or  Grx  as  given  by  equation  (60)),  which 
may  be  very  large  for  flows  generated  by  centrifugal  forces,  in  the 
denominator  of  the  ordinate  implies  that  velocities  of  appreciable  mag¬ 
nitude  can  be  associated  with  such  free -convection  flows.  In  particular, 
if  fjr  =  106  feet  per  second  squared,  which  is  a  reasonable  conservative 
figure  for  present-day  rotating  systems,  e=  0.2  {which  is  within  the 
limits  of  the  theory  presented  herein),  and  arbitrarily  X  =0.25  foot, 
then  the  maximum  velocity  attained  at  a  Prandtl  number  of  0.72  is  approx¬ 
imately  125  feet  per  second.  This  value  of  the  maximum  velocity  could, 
of  course,  be  doubled  or  even  tripled  under  the  proper  conditions.  One 
limitation  to  a  calculation  of  this  sort,  as  can  be  seen  by  comparison 
of  the  denominators  of  the  left  and  middle  terms  of  equation  (60), 
should  be  kept  in  mind;  namely,  the  limiting  Grashof  number  for  laminar 
flows.  In  lieu  of  a  complete  stability  analysis  on  this  type  of  flow, 
this  limiting  value  is  taken  to  be  10^,  as  indicated  in  reference  10. 
Consideration  of  this  limitation  then  implies  (see  equation  (60))  that 
for  large  laminar  velocities  either  must  be  large  or  X  must  be 

small. 


Comparison  with  Experiments 

Careful  experiments  of  free -convection  flows  (as  generated  by  gravi¬ 
tational  forces)  about  vertical  flat  plates  were  made  by  Schmidt  and 
Beckmann  (reference  2)  in  which  velocity  measurements  at  various  points 
along  the  plate  were  made  by  means  of  a  quartz -filament  anemometer  and 
the  temperature  measurements  were  obtained  by  means  of  manganese - 
constantan  thermocouples.  Eckert  (reference  l)  performed  similar  experi¬ 
ments  in  which  the  measurements  were  made  by  means  of  a  Zehnder-Mach 
interferometer.  The  results  of  both  sets  of  experiments  are  in  good 
agreement  with  each  other,  but  since  the  data  presented  in  reference  2 
by  Schmidt  and  Beckmann  appear  in  more  detail,  these  data  will  be  used 
for  comparison  with  the  theory. 

The  experiments  of  reference  2  were  performed  on  two  different  (in 
that  the  edges  were  smoothed  either  symmetrically  or  not)  12-  by  25-centi¬ 
meter  plates  and  on  one  50-  by  50-centimeter  plate.  It  should  here  be 
pointed  out  that  the  results  for  the  two  smaller  plates  were  almost  iden¬ 
tical  and  that  the  flow  was  entirely  laminar  except  near  the  outer  edge 
of  the  boundary  layer  where  the  slight  turbulence  of  the  room  air  dis¬ 
turbed  the  measurements  somewhat.  (This  effect  was  also  observed  by 
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Eckert.)  Large  periodic  oscillations  of  the  flow  near  the  downstream 
edge  of  the  larger  plate  were  observed  in  addition  to  the  slight  turbu¬ 
lence  near  the  outer  edge  of  the  boundary  layer.  Hence  the  data  from 
the  larger  plate  should  not  be  expected  to  yield  completely  satisfactory 
agreement  with  the  laminar  theory  as  presented  here. 

Since  the  physical  quantities  can  be  expressed  in  terms  of  a  single 
variable  as  in  equations  (60)  and  (61),  it  is  to  be  expected  that  the 
data  taken  at  the  various  points  along  the  plates  should  all  lie  on  a 
single  line  if  the  data  are  correlated  according  to  equations  (60)  and 
(61).  Thus  for  the  smaller  plates  where  (Tq-T«)  =  95.22°  B  and 
T»  =  518.68°  E,  equations  (60)  to  (62)  become 


U  -F’(tj) 

4 ,862-^/x 

(63) 

T  -  518 .68  -rrf-.  \ 

95.22  “  H(T|' 

(64) 

rj  =  88.26  — y- 

(65) 

X4 


The  velocity  and  temperature  distributions  are  so  plotted  in  figures  3  and  4, 
respectively,  as  are  the  curves  computed  theoretically  for  Pr  =  0.72.  It 
can  be  seen  that  the  agreement  is  in  general  very  good  for  small  values  of 
Tj  and  somewhat  less  satisfactory  though  still  rather  good  for  the  larger 
values  of  tj.  The  scatter  in  the  range  of  the  larger  values  of  tj  is 
believed  to  be  caused  by  the  previously  discussed  room  turbulence.  It 
should  also  be  noted  that  the  points  farthest  away  from  the  theoretical 
are  those  measured  near  the  leading  edge.  These  points  should  not,  of 
course,  be  expected  to  agree  too  well  with  the  theory  since  the  boundary- 
layer  assumptions  made  in  the  theoretical  development  imply  that  the  dis¬ 
tance  along  the  plate  is  large  as  compared  with  the  boundary -layer  thick¬ 
ness  .  Hence,  this  assumption  is  invalid  near  the  leading  edge.  Schmidt 
wnfl  Beckmann  obtained  closer  agreement  between  the  theory  and  the  experi¬ 
ments  for  the  temperature  data  and  poorer  agreement  for  the  velocity  data 
by  basing  the  kinematic  viscosity  coefficient  in  equation  (62)  on  the 
plate  temperature  rather  than  on  the  undisturbed  stream  temperature  as  was 
done  here . 

For  the  larger  plate,  (Tq-TJ  =  83.7°  E  and  T„  =  527.14°  E  so 
that  equations  (60)  to  (62)  become 


_ 2 -  =  F'  (tj) 

4.522-y/X 


0 


(66) 


n  j  C9 
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T  -  527.14  ,_x 

83.7  ”  S(t|)  '  (67) 

tj  *  83.93  -JL-  (68) 


The  yeloci ty  and  temperature  distributions  for  this  experiment  are  plotted 
in  figures  5  and  6,  respectively,  and  again  the  theoretical  curves  for 
Pr  =  0.72  are  included.  In  figure  5  it  can  be  seen  that  for  large  tj 
the  agreement  is  rather  poor,  particularly  for  the  data  for  both  small 
and  large  values  of  X.  The  poor  agreement  for  small  values  of  X  is 
again  due  to  the  theory  limitation  near  the  edge  of  the  plate  and  for 
the  large  values  of  X,  to  the  fact  that  the  flow  was  becoming  turbulent 
there . 


Flow  and  Heat -Transfer  Parameters 

In  addition  to  the  velocity  and  temperature  distributions,  it  is 
often  desirable  to  compute  other  physically  important  quantities  (such 
as  shear  stress,  drag,  heat -transfer  rate,  and  heat -transfer  coefficient) 
associated  with  the  free -convection  flow.  To  this  end,  two  parameters, 
a  flow  parameter  and  a  heat -transfer  parameter,  are  derived  in  appen¬ 
dixes  C  and  D,  respectively. 

The  flow  parameter 

- 1 -  -F”(0) 

(4  Grx3)4  <v-  ^o/*2> 

is  presented  as  a  function  of  Prandtl  number  in  figure  7 .  Thus  the 
various  flow  quantities  for  a  given  set  of  conditions  can  easily  be 
computed  by  application  of  figure  7 . 

The  local  heat -transfer  parameter 

— _  =  -H'(0) 

(Grx/4) 

as  determined  here  is  given  as  a  function  of  Prandtl  number  in  figure  8. 

A  calculation  of  the  local  Nusselt  number  from  this  equation  for 

Pr  =  0.72  and  Grx  =  109  yields  a  value  of  63.5,  which  indicates  that 

large  heat -transfer  coefficients  can  also  be  obtained  with  free -convection 
flows . 
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On  the  basis  of  a  simplified  theory  (that  is,  by  use  of  integrated 
momentum  and  energy  equations  and  assumed  velocity  and  temperature  dis¬ 
tributions),  Eckert  (see  p.  162  of  reference  7)  obtained  the  approximate 
relation 

1 

Nu  0.718(Pr)2 

1  ~  1 

(Grx/4)4  (0.952  +  Pr)4 


The  curve  representing  this  equation  is  also  presented  in  figure  8,  and 
it  closely  approximates  (to  within  about  10  percent)  the  curve  determined 
by  the  more  exact  considerations  of  this  report  over  the  entire  Prandtl 
number  range.  A  semiempirical  equation  as  given  in  reference  11  relating 
the  average  (over  the  length  X)  flusselt  number  to  the  Prandtl  and 
Grashof  numbers  which  has  been  used  in  the  heat-transfer  calculations  up 
to  the  present  is 


Nuay  =  0.548 


|j[Pr)GrJ 


1 

4 


The  constant  0.548  pertains  specifically  to  air;  for  oil  it  should  he 
0.555  (see  reference  12)  and  for  mercury,  approximately  0.33  (see  ref¬ 
erence  13).  In  order  to  obtain  local  values  from  the  average  ones  given 
by  the  last  equation,  it  is  merely  necessary  to  multiply  the  average 
values  by  0.75.  (The  determination  of  this  reduction  factor  of  0.75  is 
discussed  in  appendix  D.)  Thus  in  terms  of  the  local  quantities  the 
semiempirical  relation  becomes 


Nu 


0.411 


[(Pr)Grx] 


1 

4 


- - -  0.581  (Pr)4 

1 

(Grx/4)* 

The  curve  given  by  this  equation  is  also  presented  in  figure  8  and  the 
agreement  with  the  theoretical  curve  is  very  good  for  Prandtl  numbers 
near  unity,  not  so  good  for  large  Prandtl  numbers,  and  very  poor  for  the 
small  Prandtl  numbers .  Of  course,  changes  of  the  constants  in  the  semi¬ 
empirical  relation  as  previously  discussed  for  the  large  or  small  Prandtl 
number  cases  (oil  and  mercury,  respectively,  for  example)  would  cause  the 
semiempirical  curve  to  approximate  the  theoretical  curve  more  closely. 

The  values  of  the  heat -transfer  parameter  obtained  experimentally  for 
mercury  (Pr  =  0.03),  air  (Pr  =  0.72),  water  (Pr  =  7),  and  oil 
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(Pr  =  75.5,  115,  190,  224,  275,  318,  368,  and  442)  are  here  reduced  "by 
the  factor  0.75  from  the  average  values  reported.  The  value  for  mercury 
Is  an  average  taken  of  four  readings  from  a  curve  since  this  experiment 
•was  the  only  one  not  reported  in  tabular  form.  From  figure  8  it  can  he 
seen  that  all  of  the  experimental  values  except  those  for  the  oil  experi¬ 
ments  are  in  very  good  agreement  with  the  theoretically  computed  values. 
The  data  from  the  oil  experiments,  though  not  so  good,  show  reasonable 
agreement  (maximum  error  of  approximately  20  percent)  with  the  theoreti¬ 
cal  curve  and  good  agreement,  as  is  to  be  expected,  with  the  semi- 
empirical  curve.  The  difference  between  the  theoretical  values  and  the 
oil  experiment  results  can  possibly  be  due  to  the  fact  that  the  viscosity 
changes  in  oil  are  large  even  for  small  temperature  differences  or  due 
to  the  end  effects  in  the  measurements . 


CONCLUSIONS 

An  analysis  was  made  of  the  free -convection  flow  about  a  flat  plate 
oriented  in  a  direction  parallel  to  that  of  the  generating  body  force 
under  the  prime  assumption  that  the  relative  temperature  difference  is 
small.  It  was  found  that  the  Grashof  number  was  the  principal  factor 
determining  the  type-  of  flow  and  that  for  large  G-rashof  numbers  the  flow 
was  of  the  boundary -layer  type.  The  theoretical  development  was  then 
continued  to  consider  only  the  cases  of  large  Grashof  number  because 
these  are  of  most  importance  in  aeronautics. 

Velocity  and  temperature  profiles  for  Prandtl  numbers  of  0.01,  0.72, 
0.733,  1,  2,  10,  100,  and  1000  were  computed  on  the  basis  of  a  constant 
body  force  and  plate  temperature  and  agreement  with  experiments  where 
the  fluid  was  air  (Prandtl  number  of  0.72)  was  good.  It  was  also  demon¬ 
strated  that,  velocities  and  Nusselt  numbers  of  the  order  of  magnitude  of 
those  obtained  in  forced-convection  could  be  obtained  in  free-convection 
flows. 

A  flow  parameter  and  a  heat-transfer  parameter  which  are  functions 
of  the  Prandtl  number  alone  were  derived.  Calculations  of  the  important 
physical  quantities  such  as  shear  stress,  heat -transfer  rate,  and  the 
like  can  be  computed  from  these  parameters.  Values  of  the  heat-transfer 
parameter  obtained  from  an  approximate  theoretical  development  and  from 
experiments  compared  with  values  computed  from  the  present  development 
showed  good  agreement  over  a  wide  range  of  Prandtl  number  (0.01  to  1000 ) . 
It  is  shown  that  the  commonly  used  semiempirical  relation  for  the  heat- 
transfer  coefficient  will  yield  good  results  only  in  restricted  Prandtl 
number  ranges . 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  October  3,  1951 
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APPENDIX  A 
SYMBOLS 


The  following  notation  is  used  in  this  report: 
BjH 
c(n) ,  D1^n^ 


P 

fi 

fX 

Gr 


coefficients  in  numerical  differentiation  and  Integra¬ 
tion  formulas 


specific  heat  at  constant  pressure 
dimensionless  velocity  function 

components  of  "body  force  per  unit  mass,  1  =  1,  2,  3 
negative  of  X-component  of  body  force  per  unit  mass 


Grashof  number, 


„  2-p  i3 

P.  fj2  e 


GrX 

g 

H 

h 

K 

k 

l 


m,  n 
N 


Nu 


Nu 


av 


P 

Pr 


Grashof  number  based  on  X 

gravitational  force  per  unit  mass  (or  acceleration 
due  to  gravity) 

dimensionless  temperature  function 
heat -transfer  coefficient 
isothermal  compressibility  coefficient,  - 

heat -conductivity  coefficient 
characteristic  length 
arbitrary  exponents 

a  number,  defined  following  equation  (26) 

Nusselt  number,  hX/k 
average  Nusselt  number 
pressure 
Prandtl  number 


ft 


> 


i 
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gas  constant 
arbitrary  exponents 
absolute  temperature 
velocity  components,  i  =  1,  2,  3 
dimensionless  velocity  components,  i  ■  1,  2,  3 
dimensionless  velocity  component  in  x-direction 
dimensionless  velocity  component  in  y-direction 
Cartesian  coordinates,  i  =  1,  2,  3 
dimensionless  Cartesian  coordinates,  i  «=  1,  2,  3 
dimensionless  Cartesian  coordinate 
dimensionless  Cartesian  coordinate 

coefficient  of  volumetric  expansion,  p 

ratio  of  specific  heats 
Laplacian  operator 

relative  temperature  difference,  P(Tq-T„) 
similarity  variable 
dimensionless  temperature  function 
step  size  used  in  numerical  calculations 
absolute  viscosity 
kinematic  viscosity 
density 

dimensionless  pressure  function 
shear  stress 

dimensionless  density  function 


stream  function 
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Subscripts : 

J  Cartesian  tensor  and  summation  subscripts 

s  denotes  evaluation  at  static  conditions  (e  *  0) 

0  denotes  evaluation  at  plate  surface 

°*  denotes  evaluation  at  undisturbed  conditions 

Subscript  notation  is  used  to  denote  partial  differentiation. 

Superscripts : 

Primes  denote  ordinary  differentiation. 

Bars  (as  a  or  j)  denote  transformed  dimensionless  quantities. 
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APPENDIX  B 

NUMERICAL  SOLUTION  OF  SIMPLIFIED  BOUNDARY -VALUE  PROBLEM 

By  Lynn  U.  Albers 


The  method  is  presented  herein  hy  which  solutions  to  the  Boundary- 
value  problem 


F"  '  +  3FF'.  *  -  2F'2  +  H  =  0 

(Bl) 

H"  +  3  Pr  FH'  =0 

(B2) 

F(0)  =  F’(0)  =  0  H(Q)  =  1 

F'(-)  =  H(*»)  =  0 

were  obtained  for  the  cases  of  Pr  equal  to  0.01,  0.72,  22/30,  1,  2,  10, 
100,  and  1000.  This  discussion  will  enable  the  results  to  be  clearly 
evaluated  and  will  perhaps  serve  a&  a  guide  in  the  numerical  solution  of 
similar  problems. 

Each  of  the  cases  of  the  problem  has  a  solution  for  a  particular 
set  of  values  for  F,!(0)  .and  H*(0),  hereinafter  called  eigenvalues. 

The  basic  approach  to  the  problem  was  to  estimate  the  eigenvalues  and  to 
integrate  out  from  zero,  obtaining  functions  which  satisfied  equations  (Bl) 
and  (B2)  at  each  step.  The  Integration  was  continued  until  the  functions 
FT  and  behaved  in  a  fashion  inconsistent  with  the  boundary  values  at 
infinity;  for  example,  when  they  became  negative  or  diverged  to  infinity. 
Improved  estimates  of  the  eigenvalues  were  then  made  on  the  basis  of  the 
results  of  preceding  runs  and  the  process  was  repeated  successively  until 
a  solution  was  obtained. 

Modifications  required  to  overcome  specific  obstacles  will  be  dis¬ 
cussed  after  sufficient  details  of  the  basic  procedure  have  been  given. 

Then  an  evaluation  of  the  accuracy  of  the  numerical  results  will  be  made. 

The  integration  process  consists  of  two  parts,  a  starting  phase  and 
an  extension  phase.  The  starting  phase  begins  with  an  estimate  of  the 
eigenvalues  FT,(0)  and  Hr  (0)  and  a  decision  on  the  step  size  x  to 
be  used.  It  continues  with  an  iterative  process  of  alternately  com¬ 
puting  F*  !  *  and  H'  *  at  the  first  four  points  and  integrating  them 
by  five-point  formulas.  This  process  and  that  in  the  extension  phase 
are  so  arranged  that  the  differential  equations  are  satisfied  at  each 
integral  multiple  of  the  step  size. 
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The  extension  phase  used  preceding  data  to  integrate  step  "by  step 
"beyond  the  fourth  point.  Diagrams  of  "both  phases  will  "be  given  after  a 
few  preliminary  explanations. 

All  integration  formulas  used  are  "based  on  the  same  idea.  If  a 
function,  for  example,  F'*',  is  known  at  five  points,  there  is  a  unique 
fourth-degree  polynomial  which  agrees  with  it  at  these  five  points.  More¬ 
over,  if  the  successive  antiderivatives  (integrals)  F'1,  F' ,  and  F  of 
F1'1  are  known  at  one  point,  there  are  unique  fifth-,  sixth-,  and 
seventh -degree  polynomials  which  are  successive  antiderivatives  of  this 
fourth-degree  polynomial  and  which  agree  with  F'1,  F1,  and  F,  respec¬ 
tively,  at  the  one  point.  It  is  then  a  simple  algebra  problem  to  deduce 
from  the  values  of  F' ' '  at  five  points  and  F,  F',  and  F' '  at  a 
single  point  the  values  of  any  of  these  four  polynomials  at  any  point. 
These  results  will  approximate  the  functions  F,  F',  F'',  and  F1  '  1  to 
a  degree  dependent  on  step  size,  the  relative  positions  of  the  points  in 
question,  and  the  magnitude  of  the  fifth  derivative  of  F' ' '  in  the 
neighborhood  of  these  points. 

The  preceding  algebra  problem  can  be  presolved  in  all  situations  that 
arise  in  the  starting  and  extension  phases  of  the  present  problem  and 
specific  integration  formulas  may  be  deduced.  These  formulas  are  dis¬ 
cussed  in  the  next  paragraph. 


Let  F' '  '  be  denoted  at  five  successive  points  by  F^' 1 ,  F^1',  Fi,'1, 
F^' 1 ,  and  F^’’.  Ih  the  starting  phase,  these  points  are  0,  x,  2x,  3x, 
and  4x,  and  Fq,  F^,  and  F^'  are  also  known.  Then  the  five  sets  of 
formulas  required  in  the  starting  phase  are 


i  =  1,  2,  3,  4  (B3) 


i  =  1,  2,  3,  4  (B4) 


i  =  1,  2,  3,  4  (B5) 


« 


t 
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Hi  =  H0  +  ixH’  + 


0  -nT^T 


m  Z/Aij(2)HV  1  =  ^  2,  3>  4  (B6) 


'i  <3=0 


VsFo+  1 


xfo  +^fo  +4iEAij(3)Fr 

Di  <3=0 


i  k  1,  2^  3,  4 

where  the  superscripts  on  the  Aj. and  D^11^  refer  to  the  order  of 
integration. 


The  constants  A. 
tables : '  J 


and  Di^n)  may  he  read  from  the  following 


For  Aj/'1')  and  : 


0 

i 

2 

3 

4 

251 

646 

-264 

106 

-19  720 

29 

124 

24 

4 

-1  90 

27 

102 

72 

42 

-3  80 

14 

64 

24 

64 

14  45 

For  A^^  and  D.^2): 


A.  .(2) 

Aij 


\<J 

i\ 

0 

1 

2 

3 

4 

1 

367 

540 

-282 

116 

-21 

1440 

2 

53 

144 

-30 

16 

-3 

90 

3 

441 

L404 

162 

180 

-27 

480 

4 

192 

64 

0 

45 
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For  A±J ^  and  D^3); 


A  (3) 

AiJ 

»i(S) 

\  J 

i\ 

0 

1 

2 

3 

4 

1 

1017 

1070 

-618 

258 

-47 

10080 

2 

331 

664 

-240 

104 

-19 

630 

3 

1431 

3726 

-486 

450 

-81 

1120 

4 

744 

_ 

2176 

96 

384 

-40 

315 

I 


It  is  now  possible  to  diagram  the  steps  of  the  starting  phase  of 
the  integration.  If  each  har  above  a  function  denotes  an  improved  esti¬ 
mate  of  it,  and  the  first  estimates  of  Fi 1 1 ,  F* 1 1 ,  F* ' 1 ,  and  F!ir  are 

JL  C*  O 

all  equal  to  F'rf,  and  similarly  for  the  Hrf;  then  the  starting  phase 
diagrams  are  u 


(i). 


F'  ' 

0  * 


pi  1 1 

0 


pi i I  pi i ,  pi i t 

1  ’  2  ’  3  ’ 


F»  ' ' )  -*•  F  ,  F' ,  F'  ' 
4  1’  l’  1 


(This  diagram  means  that  the  values  in  parentheses  are  used  with  approp¬ 
riate  integration  formulas  from  (B3)  to  (B7)  to  obtain  F,  F',  and  F1  ' 
at  t)  =  x. ) 

(2).  (H0,  Hq,  Hq'  ,  H”,  H”,  H”, 


(3). 


F>', 


H„  .  H 


PI  I  t 
F1  ' 


(The  preceding  diagram  means  that  the  values  in  parentheses  are  substi¬ 
tuted  in  the  differential  equations  (Bl)  and  (B2)  to  obtain  F'  *  !  and 


1'  » 

at  rj 

=  X 

.) 

(4). 

<v 

Fi> 

Fo’> 

Fo' 

'  .  F' 

;  2 

II  TJM 

9  F3 

ii  Tp  i  t  t  ^  F 

9  *4  ’  29 

F’,  F- 

(5). 

<H0, 

Eo> 

E0> 

si’ 

>E2> 

TT'  ' 
*3  ’ 

Si'>-  V  H2 

(6). 

(F2’ 

F2> 

F2 '  ’ 

E2’ 

H-)- 

pi  i  t 

2 

H1  1 
>  n2 

(7). 

(F0> 

F6’ 

F6’> 

Fo 

F-' 

i  pi 

9  2 

II  Tp  1  ' 

",  F'")  -  F3, 

F3j  F3 

it 
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A 


* 


(8). 

(9) . 

(10) 
(11) 
(12) 

“been 
1  to 

(13) 

(14) 

(15) 

(16) 


(=0. 

H6, 

E6’' 

sr 

>  E2^ 

E3 '  * 

=;■')  - 

E3*  =3  * 

<F3’ 

F3' 

f>-, 

H3> 

EP  - 

F3 

H'  * 

»  3 

(F0, 

Fi7 

F0’ 

Fo' 

ri  pi  i 

7  3 

•,  F-")  -F4,  F«,  F- 

<K0> 

H6> 

H'  ' 

o  * 

Ei’ 

H T  * 

>  *2  > 

H*  1 
o  > 

=4’)^ 

H4;  E4  ' 

F4> 

E4> 

^)-* 

T'  1  1 

4 

H'  ' 

7  n4 

It  may  be  noted  'here  that  all  four  values  of  F1 1 '  and  H’ 1  have 
improved,  and  further  improvement  will  require  iteration  of  steps 
12.  The  start  of  the  second  iteration  is  diagrammed  as  follows; 


(F  , 

o' 

F-', 

F1",  F'  "  ,  F' 
0  7  1  7  2 

",  f"  '  ,  F'  "  )  -*• 

7  3  4 

V 

Fi7 

jit  ! 
1 

(H0, 

H6> 

E6"7 

5^,  H-,  H-’, 

E4'^EP  Ei 

<?1' 

Fi> 

Fj_’, 

Hi,  E-,)- 

f- 

(F07 

Fo> 

■pt  ! 

0  ' 

Fq '  1 ,  F’ 

",  F£",  !”")-»• 

F2> 

Fr  1 

2 

• 

• 

• 

• 

Successive  sets  of  12  steps  are  performed  until  the  values  of  Fj' ' 
and  HI '  no  longer  change. 

On  the  IBM  Card  Programmed  Electronic  Calculator,  a  deck  of  punched 
cards  2  inches  thick  sufficed  to  perform  steps  1  to  12 .  Three  runs  of 
this  starter  deck  at  3  minutes  per  run  accomplished  complete  convergence 
in  most  cases.  At  the  end  of  the  starting  process  there  have  teen  com¬ 
puted  and  stored  F4,  F£,  F “ ,  H4,  and  H£,  and  final  estimates  of  Fj", 

F£",  F^",  F^",  H”,  H”,  H£>,  and  H^' . 

The  extension  phase  has  now  teen  reached.  It  used  a  different  set 
of  integration  formulas  tased  on  the  same  general  ideas  as  equations  (B3) 
to  (B7)  .  If  F^",  F[",  F£",  F^",  and  '  now  designate  F'  1 '  at 

any  five  successive  points,  and  the  subscript  5  denotes  the  next  point, 
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(B8) 


(B9) 


(BIO) 


4 

%-H4  +  >HA  +  T5xEBJ<2>Hj'  (B11) 

c  j=0 

4 

F5  =  F4  +  XF|  +  it  FI'  +  ^  b/3)F'  '  ’  (B12) 

C  j=0 

•where  the  B.^n^  and  are  given  in  the  following  table: 


Bi(n) 

C(n) 

\ 

n 

j 

0 

1 

2 

3 

4 

1 

251 

-1274 

2616 

-2774 

1901 

720 

2 

135 

-692 

1446 

-1596 

1427 

1440 

3 

410 

-2116 

4476 

-5084 

5674 

20160 

The  extension  phase  may  then  "be  diagrammed  simply  as  follows: 

(1)  (f4,  F',  F-',  F-',  F'",  F-'S  F>",  f”’)-  f5,  f«,  F-' 

(2)  (H4,  H-,  H”,  I",  H”,  H”,  H-')-  Hg,  H' 


(3)  (F5,  F,L,  F£*,  H5,  H’)-*  f>",  H 


A 


A 
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The  values  of  the  functions  at  the  next  point  are  computed  in  simi¬ 
lar  manner ,  where  the  latest  sets  of  five  values  of  F1  1  1  and  Hrf  are 
used.  This  process  advances  step  by  step  toward  infinity. 

The  extender  deck  of  punched  cards  was  about  3  inches  thick,  and 
took  a  little  over  3  minutes  per  run.  For  Pr  =  0.72,  a  step  size  of 
0.1  was  used,  the  starting  phase  took  10  minutes,  and  the  extension  phase 
about  30  minutes.  When  it  is  realized  that  about  11,000  operations  were 
performed  in  the  40  minutes  per  run,  it  may  be  seen  that  solution  of  the 
present  problem  would  have  been  prohibitively  difficult  on  desk-type  cal¬ 
culators.  Simplifications  in  method  would  have  sacrificed  accuracy  or 
required  smaller  step  size. 

In  two-point  boundary  value  problems  where  one  point  is  infinity, 
some  problems  of  judgment  are  involved  as  to  where  infinity  is,  and  as  to 
when  a  satisfactory  approximation  to  a  solution  has  been  obtained.  In 
most  cases  this  question  was  settled  for  the  present  problem  by  calling 
a  run  satisfactory  when  it  fell  between  two  runs  for  which  FT  and  H 
did  not  differ  at  important  points  in  the  fourth  decimal  place,  and  for 
which  F1  and  H  flattened  out  at  zero,  correct  to  four  decimal  places. 


Certain  difficulties  were  met  in  the  attempt  to  use  the  basic  pro¬ 
cedure  previously  discussed.  These  necessitated  certain  modifications. 

For  Pr  =  2,  10,  100,  and  1000,  H  would  settle  down  to  zero  at  an 
early  stage,  but  while  F*  was  still  coming  down,  H,!  would  begin  to 
oscillate  and  these  oscillations  increased  and  fed  back  into  all  other 
functions.  This  trouble  was  avoided  by  the  following  modifications:  It 
is  a  consequence  of  equation  (B2)  that 


H'  (t))  =  H'  (0)exp  ^3  Pr  J"  F(t)  dtj 


=  H'  (t)  -x  )exp 


^■3  Pr F(t)  dt^ 


(B13) 


The  extension  phase  was  modified  to  require  the  additional  integra¬ 
tion  formulas 


E 

i=l 


Ai  Fi 


(B14) 
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Hi 


A 

(B15) 


where  Ax  =  -19,  Ag  =  106,  A3  =  -264,  A 4  =  646,  and  Ag  =  251. 

These  formulas  were  used  along  with  equations  (B8)  to  (BIO)  according 
to  the  following  diagram: 


(D  (F4,  F',  F”,  F'”,  F*",  F-",  F-",  F--«)  -  Fg,  F',  F” 

(2)  (F1,  F2,  F3,  F^,  F5,  Hp  -+  Hg  hy  means  of  (B14)  and  (B13) 


(3)  (Hp  H£,  H|,  H£,  H4)  ■+  H5  hy  means  of  (B15) 

(4)  (F5,  F,L,  Fg* ,  H5)  -  F^  * 


The  value  F^tT  is  discarded  and  Fr  1 1  at  the  last  five  points  is 

used  to  repeat  the  whole  process  again  and  again  ad  infinitum*  As  long 
as  F  stays  positive,  HT  is  guaranteed  to  approach  zero  and  H  will 
flatten  out  to  some  value  and  not  oscillate. 


For  Fr  =  0.01,  0.72,  22/30,  and  1,  the  Fr  T  T  "began  to  oscillate 
at  an  advanced  point  and  these  oscillations  grew  and  fed  into  the  other 
functions.  For  all  cases  hut  Pr  =  0.01,  the  oscillations  appeared  very 
late,  hear  the  end  of  the  run,  and  a  suitable  halving  of  step  size  when 
oscillation  was  detected  in  the  fourth  differences  of  FItf  was  suffic¬ 
ient  to  avoid  the  difficulty.  But  for  the  0.01  case,  oscillations  of 
FT,T  appeared  early  in  the  run,  namely,  soon  after  the  peak  in  F. 

These  oscillations  were  found  to  he  step-size  connected,  so  that  reduc¬ 
tion  of  the  step  to  0.02  avoided  them.  Even  then  oscillations  in  F,,f 
would  begin  to  appear  e^ery  25  steps  or  so,  and  these  were  smoothed  out 
regularly  hy  repeated  runs  of  a  deck  similar  to  the  starter  deck.  Each 
run  under  these  conditions  took  about  16  hours,  making  this  the  most 
difficult  case  to  solve. 
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APPENDIX  C 

DERIVATION  OF  FLOW  PARAMETER 


By  definition  the  shear  stress  is  given  hy 

faj\ 

T ' 


(01) 


To  express  (Su/By)q  in  terms  of  the  known  function  F(r|), 
made  of  equations  (60)  and  (62).  Then 


1 

4  v„ 


use  can  "be 


Substitution  of  this  expression  into  equation  (Cl)  yields  the  flow 
parameter 


- - =  F"(0) 

(4  ,r/)\  ^?) 

Note  that  from  the  general  derivation,  the  flow  parameter  contains  the 
.viscosity  evaluated  at  two  different  points.  Recall,  however,  that  the 
analysis  has  shown  that  to  a  first  approximation  the  variation  of  vis¬ 
cosity  with  temperature  can  he  neglected.  Thus  the  viscosity  can  he 
taken  as  constant  in  the  entire  flow  field. 
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appendix:  d 

DERIVATION  OP  EEAT -TRANSFER  PARAMETER 


The  local  Nusselt  number  is  defined  aa 


„T  _  hX  -X  /dT\ 

k  (Tq-^)  \dj)Q 


(Dl) 


To  express  (St/By)q 
of  equations  (61)  and 


in  terms  of  the  known  function 
(62) .  Thus 


to  <To-E.)  (GrxY  , 

S' — s— w 


H(n)  ,  use  is  made 


Substitution  of  this  expression  into  equation  (Dl)  yields  the  heat- 
transfer  parameter 


A 


- — — j-  =  -  H'  (0)  (D2) 

(Grx/4)4 

The  heat -transfer  parameter  as  given  by  equation  (D2)  is,  as  was  pre¬ 
viously  stated,  a  local  parameter.  It  is  often  desired  to  compute  the 
average  (over  the  length  X)  value  of  this  parameter.  To  this  end,  the 
Nusselt  number  (as  given  in  equation  (Dl))  must  be  defined  in  terms  of 
an  average  heat-transfer  coefficient  and  the  quantity  thus  obtained  must 
then  be  integrated  over  the  length  X  and  divided  by  X.  This  procedure 
yields  the  result 

Nu  =  t-  (Nu) 

4  N  'av 

It  is  from  this  last  equation  that  the  0.75  reduction  factor  previously 
discussed  was  obtained. 
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TABLE  I  -  FUNCTIONS  F  AND  H  AND  DERIVATIVES  FOE  VARIOUS  PRANDTL  NUMBERS 

(a)  Prandtl  number,  O.OI 


n 

F 

F* 

Fr ' 

E 

E* 

n 

F 

F- 

F‘  » 

H 

H' 

0 

0.0000 

0.0000 

0.9862 

1.0000 

-0.0812 

4.16 

1.9319 

0.4037 

-0.0615 

0.6741 

-0.0723 

.i 

.0048 

.0936 

.8868 

.9919 

-.0812 

4.18 

1.9400 

.4024 

-.0614 

.6726 

-.0722 

.2 

.0184 

.1774 

.7891 

.9838 

-.0812 

4.20 

1.9480 

.4012 

-.0612 

.6712 

-.0721 

.3 

.0399 

.2516 

.6942 

.9756 

-.0812 

4.24 

1.9640 

.3988 

-.0609 

.6683 

-.0720 

.4 

.0684 

.3164 

.6030 

.9675 

-.0811 

4.28 

1.9799 

.3964 

-.0605 

.6654 

-.0718 

.5 

.1029 

.3723 

.5163 

.9594 

-.0811 

4.32 

1.9957 

.3939 

-.0602 

.6626 

-.0716 

.6 

.1426 

.4198 

.4349 

.9513 

-.0811 

4.36 

2.0114 

.3915 

-.0599 

.6597 

-.0714 

.7 

.1866 

.4595 

.3595 

.9432 

-.0811 

4.40 

2.0270 

.3891 

-.0595 

.6568 

-.0713 

.8 

.2342 

.4919 

.2906 

.9351 

-.0810 

4.50 

2.0657 

.3832 

-.0587 

.6497 

-.0708 

.9 

.2848 

.5178 

.2285 

.9270 

-.0809 

4.60 

2.1037 

.3774 

-.0579 

.6427 

-.0704 

1.0 

.3376 

.5379 

.1734 

.9189 

-.0809 

4.70 

2.1411 

.3716 

-.0572 

.6357 

-.0699 

1.1 

.3922 

.5527 

.1253 

.9108 

-.0808 

4.80 

2.1780 

.3660 

-.0564 

.6287 

-.0695 

1.2 

.4480 

.5631 

.0839 

.9028 

-.0807 

4.90 

2.2143 

.3604 

-.0556 

.6218 

-.0690 

1.3 

.5047 

.5697 

.0489 

.8947 

-.0806 

5.00 

2.2501 

.3548 

-.0549 

.6149 

-.0686 

1.4 

.5618 

.5731 

.0198 

.8866 

-.0804 

5.10 

2.2853 

.3494 

-.0541 

.6081 

-.0681 

1.5 

.6192 

.5739 

-.0039 

.8786 

-.0803 

5.20 

2.3200 

.3440 

-.0534 

.6013 

-.0676 

1.6 

.6765 

.5725 

-.0229 

.8706 

-.0801 

5.40 

2.3877 

.3335 

-.0520 

.5878 

-.0667 

1.7 

.7336 

.5694 

-.0379 

.8626 

-.0800 

5.60 

2.4534 

.3232 

-.0506 

.5746 

-.0657 

1.8 

.7904 

.5650 

-.0493 

.8546 

-.0798 

5.80 

2.5170 

.3132 

-.0493 

.5615 

-.0648 

1.9 

.8466 

.5596 

-.0579 

.8466 

-.0796 

6.00 

2.5787 

.3035 

-.0479 

.5487 

-.0638 

2.0 

.9023 

.5535 

-.0641 

.8387 

-.0794 

6.20 

2.6384 

.2940 

-.0467 

.5360 

-.0628 

2.1 

.9573 

.5469 

-.0685 

.8307 

-.0792 

6.40 

2.6963 

.2849 

-.0454 

.5236 

-.0618 

2.2 

1.0117 

.5399 

-.0714 

.8228 

-.0789 

6.60 

2.7524 

.2759 

-.0442 

.5113 

-.0608 

2.3 

1.0653 

.5326 

-.0732 

.8150 

-.0787 

6.80 

2.8067 

.2671 

-.0430 

.4993 

-.0598 

2.4 

1.1182 

.5253 

-.0742 

.8071 

-.0784 

7.00 

2.8593 

.2586 

-.0419 

.4874 

-.0588 

2.5 

1.1703 

.5178 

-.0746 

.7993 

-.0782 

7.40 

2.9594 

.2423 

-.0397 

.4643 

-.0568 

2.6 

1.2217 

.5104 

-.0745 

.7915 

-.0779 

7.80 

3.0532 

.2269 

-.0376 

.4420 

-.0547 

2.7 

1.2724 

.5029 

-.0741 

.7837 

-.0776 

8.20 

3.1411 

.2123 

-.0356 

.4205 

-.0527 

2.8 

1.3223 

.4956 

-.0735 

.7760 

-.0773 

8.60 

3.2232 

.1984 

-.0337 

.3998 

-.0508 

2.9 

1.3715 

.4882 

-.0728 

.7682 

-.0770 

9.00 

3.2999 

.1853 

-.0319 

.3799 

-.0489 

3.0 

1.4200 

.4810 

-.0719 

.7606 

-.0766 

9.40 

3.3715 

.1729 

-.0302 

.3607 

-.0470 

3.1 

1.4677 

.4739 

-.0711 

.7529 

-.0763 

9.80 

3.4383 

.1612 

-.0285 

.3423 

-.0451 

3.2 

1.5148 

.4668 

-.0701 

.7453 

-.0760 

10.20 

3.5005 

.1501 

-.0270 

.3246 

-.0432 

3.3 

1.5611 

.4598 

-.0692 

.7377 

-.0756 

10.60 

3.5584 

.1396 

-.0254 

.3069 

-.0413 

3.4 

1.6067 

.4530 

-.0683 

.7302 

-.0753 

11.00 

3.6123 

.1297 

-.0241 

.2915 

-.0397 

3.45 

1.6293 

.4496 

-.0678 

.7264 

-.0751 

11.40 

3.6622 

.1203 

-.0228 

.2759 

-.0380 

3.5 

1.6517 

.4462 

-.0674 

.7227 

-.0749 

11.80 

3.7086 

.1114 

-.0215 

.2610 

-.0364 

3.55 

1.6739 

.4428 

-.0669 

.7189 

-.0747 

12.20 

3.7514 

.1031 

-.0203 

.2468 

-.0348 

3.6 

1.6960 

.4395 

-.0664 

.7152 

-.0745 

12.60 

3.7911 

.0952 

-.0192 

.2332 

-.0332 

3.65 

1.7179 

.4362 

-.0660 

.7115 

-.0743 

13.00 

3.8276 

.0877 

-.0181 

.2202 

-.0317 

3.7 

1.7396 

.4329 

-.0655 

.7078 

-.0741 

13.60 

3.8771 

.0773 

-.0165 

.2018 

-.0296 

3.75 

1.7611 

.4296 

-.0651 

.7041 

-.0739 

14.20 

3.9206 

.0679 

-.0151 

.1847 

-.0276 

3.8 

1.7825 

.4264 

-.0646 

.7004 

-.0738 

14.80 

3.9587 

.0592 

-.0138 

.1687 

-.0257 

3.85 

1.8308 

.4232 

-.0642 

.6967 

-.0736 

15.40 

3.9918 

.0513 

-.0126 

.1538 

-.0239 

3.9 

1.8249 

.4200 

-.0638 

.6930 

-.0734 

16.00 

4.0204 

.0441 

-.0114 

.1399 

-.0223 

3.95 

1.8458 

.4168 

-.0633 

.6893 

-.0732 

17.00 

4.0591 

.03  3  5 

-.0097 

.1189 

-.0197 

4.00 

1.8665 

.4136 

-.0629 

.6857 

-.0729 

18.00 

4.0880 

.0246 

-.0082 

.1003 

-.0175 

4.05 

1.8871 

.4105 

-.0625 

.6821 

-.0727 

19.00 

4.1088 

.0171 

-.0069 

.0839 

-.0154 

4.10 

1.9076 

.4074 

-.0620 

.6784 

-.0725 

20.00 

4.1226 

.0108 

-.0057 

.0694 

-.0137 

4.12 

1.9157 

.4061 

-.0619 

.6770 

-.0725 

21.00 

4.1308 

.0057 

-.0046 

.0565 

-.0121 

4.14 

1.9238 

.4049 

-.0617 

.6755 

-.0724 

22.00 

4.1343 

.0015 

-.0037 

.0452 

-.0107 
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(t)  Prandtl  number,  0.72 

tj  F  F’  ]  FT  1  H  |  H' 

0  0.0000  0.0000  0.6760  1.0000  -0.5046 

.1  .0032  .0627  .5785  .9495  -.5045 

.2  .0122  .1159  .4866  .8991  -.5037 

.3  .0261  .1602  .4007  .8488  -.5016 

.4  .0440  .1962  .3210  .7989  -.4979 

.5  .0651  .2246  .2479  .7493  -.4921 

.6  .0887  .2460  .1817  .7005  -.4840 

.7  .1141  .2612  .1224  .6526  -.4735 

.1407  .2708  .0700  .6059  -.4607 

.1681  .2754  .0246  .5606  -.4456 

.1957  .2759  -.0143  .5168  -.4284 

.2231  .2728  -.0468  .4749  -.4095 

1.2  .2501  .2667  -.0734  .4350  -.3891 

1.3  .2764  .2583  -.0945  .3972  -.3676 

1.4  .3017  .2480  -.1106  .3615  -.3453 

1.5  .3260  .2363  -.1224  .3281  -.3227 

1.6  .3490  .2236  -.1302  .2970  -.3000 

1.7  .3707  .2104  -.1347  .2681  -.2775 

1.8  .3910  .1968  -.1363  .2415  -.2556 

1.9  .4100  .1832  -.1356  .2170  -.2344 

2.0  .4277  .1697  -.1331  .1945  -.2141 

2.1  .4440  .1566  -.1290  .1741  -.1949 

2.2  .4590  .1440  -.1238  .1555  -.1768 

2.3  .4728  .1319  -.1178  .1387  -.1598 

2.4  .4854  .1204  -.1113  .1235  -.1441 

2.5  .4969  .1097  -.1044  .1098  -.1296 

.2.6  .5074  .0996  -.0973  .0976  -.1163 

2.7  .5168  .0902  -.0903  .0865  -.1041 

2.8  .5254  .0815  -.0834  .0767  -.0930 

2.9  .5332  .0735  -.0767  .0679  -.0830 

3.0  .5401  .0661  -.0703  .0601  -.0739 

3.1  .5464  .0594  -.0642  .0531  -.0657 

3.2  .5520  .0533  -.0585  .0469  -.0584 

3.3  .5571  .0477  -.0531  .0414  -.0518 

3.4  .5616  .0427  -.0481  .0365  -.0459 

3.6  .5692  .0339  -.0392  .0284  -.0359 

3.8  .5735  .0269  -.0317  .0220  -.0281 

4.0  .5801  .0212  -.0255  .0170  -.0219 

4.2  .5838  .0166  -.0204  .0132  -.0170 

4.4  .5868  .0130  -.0162  .0102  -.0132 

4.6  .5891  .0101  -.0128  .0078  -.0102 

4.8  .5908  .0078  -.0102  .0060  -.0080 

5.0  .5922  .0060  -.0079  .0046  -.0061 

5.2  .5933  .0046  -.0063  .0035  -.0048 

5.4  ,5941  .0034  -.0048  .0027  -.0037 

5.8  .5951  .0019  -.0029  .0015  -.0022 

6.2  .5957  .0010  -.0018  .0009  -.0013 

6.8  .5960  .0003  -.0008  .0003  -.0006 

7.3  .5961  .0000  -.0004  .0001  -.0003 


TABLE  I  -  FUNCTIONS  F  AND  H  AND  DERIVATIVES 


(c)  Prandtl  number,  0.733 


n 

F 

F  ’  F1  ' 

H 

H* 

0 

0.0000 

0.0000  0.6741 

1.0000 

-0.5080 

.i 

.0032 

.0625  .5767 

.9492 

-.5079 

.2 

.0122 

.1155  .4849 

.8984 

-.5071 

.3 

.0260 

.1597  .3990 

.8478 

-.5050 

.4 

.0438 

.1955  .3194 

.7975 

-.5012 

.5 

.•0649 

.2238  .2465 

.7477 

-.4953 

.6 

.0884 

.2451  .1804 

.6985 

-.4870 

.7 

.1137 

.2601  .1212 

.6503 

-.4763 

.8 

.1402 

.2695  .0691 

.6033 

-.4632 

.9 

.1674 

.2741  .0237 

.5578 

-.4478 

1.0 

.1949 

.2745  -.0149 

.5139 

-.4303 

1.1 

.2222 

.2713  -.0473 

.4718 

-.4110 

1.2 

.2490 

.2652  -.0737 

.4317 

-.3902 

1.3 

.2752 

.2568  -.0946 

.3938 

-.3684 

1.4 

.3003 

.2465  -.1106 

.3581 

-.3458 

1.5 

.3244 

.2348  -.1222 

.3246 

-.3228 

1.6 

.3473 

.2222  -.1299 

.2935 

-.2998 

1.7 

.3688 

.2090  -.1342 

.2647 

-.2771 

1.8 

.3891 

.1954  -.1358 

.2381 

-.2549 

1.9 

.4079 

.1819  -.1350 

.2136 

-.2335 

2.0 

.4254 

.1685  -.1324 

.1913 

-.2130 

2.1 

.4416 

.1554  -.1283 

.1710 

-.1937 

2.2 

.4565 

.1429  -.1230 

.1526 

-.1754 

2.3 

.4702 

.1309  -.1170 

.1359 

-.1584 

2.4 

.4827 

.1195  -.1104 

.1208 

-.1427 

2.5 

.4941 

.1088  -.1035 

.1073 

-.1281 

2.6 

.5045 

.0988  -.0965 

.0952 

-.1148 

2.7 

.5139 

.0895  -.0895 

.0843 

-.1026 

2.8 

.5224 

.0809  -.0826 

.0746 

-.0916 

2.9 

.5301 

.0729  -.0759 

.0660 

-.0816 

3.0 

.5370 

.0657  -.0695 

.0583 

-.0725 

3.1 

.5433 

.0590  -.0635 

.0514 

-.0644 

3.2 

.5489 

.0530  -.0578 

.0454 

-.0571 

3.3 

.5539 

.0475  -.0524 

.0400 

-.0506 

3.4 

.5584 

.0425  -.0474 

.0352 

-.0448 

3.6 

.5660 

.0339  -.0386 

.0273 

-.0350 

3.8 

.5720 

.0269  -.0312 

.0211 

-.0272 

4.0 

.5769 

.0213  -.0250 

.0163 

-.0211 

4.2 

.5807 

.0169  -.0200 

.0126 

-.0164 

4.4 

.5837 

.0133  -.0159 

.0097 

-.0127 

4.6 

.5860 

.0105  -.0125 

.0074 

-.0098 

4.8 

.5879 

.0082  -.0099 

.0057 

-.0076 

5.0 

.5893 

.0065  -.0077 

.0044 

-.0058 

5.2 

.5905 

.0051  -.0061 

.0033 

-.0045 

5.4 

.5914 

.0040  -.0047 

.0025 

-.0035 

5.6 

.5921 

.0032  -.0036 

.0019 

-.0027 

5.8 

.5927 

.0025  -.0028 

.0015 

-.0021 

6.0 

.5932 

.0021  -.0022 

.0011 

-.0016 

6.4 

.5938 

.0014  -.0013 

.0006 

-.0009 

6.8 

.5943 

.0010  -.0008 

.0003 

-.0006 

7.2 

.5946 

.0007  -.0004 

.0001 

-.0003 

7.6 

.5949 

.0006  -.0002 

.0000 

-.0002 

8.0 

.5951 

.0005  -.0001 

.0000 

-.0001 
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FOR  VARIOUS  RRAKDTL  NUMBERS  -  CONTINUED 


(d)  Prandtl  number,  1 


n 

F 

pt 

F'  ' 

H 

H' 

0 

0.0000 

o 

o 

o 

o 

o 

0.6421 

1.0000 

-0.5671 

.i 

.0030 

.0593 

.5450 

.9433 

-.5669 

.2 

.0115 

.1092 

.4540 

.8867 

- .5657 

.3 

.0246 

.1503 

.3694 

.8302 

-.5627 

.4 

.0413 

.1833 

.2916 

.7742 

-.5572 

.5 

.0610 

.2089 

.2208 

.7189 

-.5488 

.6 

.0829 

.2277 

.1572 

.6645 

-.5371 

.7 

.1064 

.2406 

.1008 

.6116 

-.5221 

.8 

.1308 

.2481 

.0516 

.5602 

-.5038 

.9 

.1558 

.2511 

.0093 

.5109, 

-.4826 

1.0 

.1809 

.2502 

-.0263 

.4638 

-.4589 

1.1 

.2058 

.2461 

-.0557 

.4192 

-.4330 

1.2 

.2300 

.2393 

-.0793 

.3772 

-.4056 

.1.3 

.2535 

.2304 

-.0975 

.3381 

-.3772 

1.4 

.2761 

.2199 

-.1110 

.3018 

-.3484 

1.5 

.2975 

.2083 

-.1203 

.2684 

-.3197  - 

1.6 

.3177 

.1960 

-.1260 

.2379 

-.2915 

1.7 

.3367 

.1832 

-.1287 

.2101 

-.2642 

1.8 

• .3543 

.1703 

-.1288 

.1850 

-.2382 

1.9 

.3707 

.1575 

-.1268 

.1624 

-.2136 

2.0 

.3859 

.1450 

-.1233 

.1422 

-.1907 

2.1 

.3997 

.1329 

-.1185 

.1242 

-.1695 

2.2 

.4125 

.1213 

-.1127 

.1082 

-.1501 

2.3 

.4240 

.1104 

-.1064 

.0941 

-.1324 

2.4 

.4346 

.1001 

-.0997 

.0817 

-.1164 

2.5 

.4441 

.0904 

-.0928 

.0708 

-.1020 

2.6 

.4527 

.0815 

-.0859 

.0613 

-.0892 

2.7 

.4604 

.0733 

-.0791 

.0529 

-.0777 

2.8 

.4673 

.0657 

-.0725 

.0457 

-.0676 

2.9 

.4736 

.0588 

-.0662 

.0392 

-.0587 

3.0 

.4791 

.0524 

-.0602 

.0339 

-.0509 

3.1 

.4841 

.0467 

-.0546 

.0291 

-.0441 

3.2 

.4885 

.0415 

-.0493 

.0250 

-.0381 

3.3 

.4924 

.0368 

-.0444 

.0215 

-.0329 

3.4 

.4959 

.0326 

-.0399 

.0185 

-.0283 

3.6 

.5016 

.0254 

-.0321 

.0136 

-.0210 

3.8 

.5061 

.0197 

-.0255 

.0099 

-.0155 

4.0 

.5096 

.0151 

-.0202 

.0072 

-.0115 

4.2 

.5122 

.0116 

-.0158 

.0053 

-.0084 

4.4 

.5143 

.0087 

-.0124 

.0038 

-.0062 

4.6 

.5158 

.0066 

-.0096 

.0027 

-.0045 

4.8 

.5169 

.0049 

-.0075 

.0020 

-.0034 

5.0 

.5177 

.0035 

-.0057 

.0014 

-.0024 

5.2 

.5183 

.0025 

-.0044 

.0010 

-.0018 

5.5 

.5189 

.0014 

-.0029 

.0006 

-.0011 

6.0 

.5194 

.0004 

-.0014 

.0002 

-.0005 

6.25 

.5194 

.0000 

-.0010 

.0000 

-.0004 

a 


(e)  Pran&tl  number,  2 


n 

F 

F’ 

pi  i 

H 

H' 

0 

o 

o 

o 

o 

o 

0 .0000 

0.5713 

1.0000 

-0.7165 

.1 

.0027 

.0523 

.4749 

.9284 

-.7161 

.2 

.0101 

.0952 

.3861 

.8569 

-.7135 

.3 

.0215 

.1297 

.3049 

.7858 

-.7069 

.4 

.0358 

.1565 

.2318 

.7157 

-.6949 

.5 

.0525 

.1763 

.1666 

.6470 

-.6768 

.6 

.0709 

.1901 

.1095 

.5805 

-.6523 

.7 

.0904 

.1985 

.0602 

.5168 

-.6215 

.8 

.1104 

,2024 

.0185 

.4564 

-.5852 

.9 

.1307 

.2024 

-.0161 

.3999 

-.5443 

1.0 

.1508 

.1994 

-.0440 

.3476 

-.5002 

1.1 

.1705 

.1938 

-.0659 

.2999 

-.4543 

1.2 

.1895 

.1864 

-.0824 

.2568 

-.4077 

1.3 

.2077 

.1775 

-.0942 

.2183 

-.3619 

1.4 

.2250 

.1677 

-.1020 

.1844 

-.3178 

1.5 

.2412 

.1572 

-.1063 

.1547 

-.2763 

1.6 

.2564 

.1465 

-.1078 

.1290 

-.2380 

1.7 

.2706 

.1357 

-.1071 

.1069 

-.2032 

1.8 

.2836 

.1252 

-.1046 

.0882 

-.1721 

1.9 

.2956 

.1149 

-.1008 

.0724 

-.1446 

2.0 

.3066 

.1050 

-.0960 

.0592 

-.1207 

2.1 

.3166 

.0957 

-.0906 

.0482 

-.1001 

2.2 

.3257 

.0869 

-.0849 

.0391 

-.0826 

2.3 

.3340 

.0787 

-.0789 

.0316 

-.0677 

2.4 

.3415 

.0711 

-.0729 

.0254 

-.0553 

2.5 

.3483 

.0641 

-.0671 

.0204 

-.0450 

2.6 

.3543 

.0577 

-.0614 

.0164 

-.0364 

2.7 

.3598 

.0518 

-.0560 

.0131 

-.0294 

2.8 

.3647 

.0465 

-.0509 

.0105 

-.0237 

2.9 

.3691 

.0417 

-.0461 

.0083 

-.0190 

3.0 

.3731 

.0373 

-.0416 

.0066 

-.0152 

3,1 

.3766 

.0333 

-.0375 

.0053 

-.0121 

3.2 

.3798 

.0298 

-.0338 

.0042 

-.0097 

3.3 

.3826 

.0266 

-.0303 

.0033 

-.0077 

3.4 

.3851 

.0237 

-.0272 

.0026 

-.0061 

3.6 

.3893 

.0188 

-.0218 

.0017 

-.0038 

3.8 

.3927 

.0149 

-.0174 

.0010 

-.0024 

4.0 

.3953 

.0118 

-.0138 

.0007 

-.0015 

4.2 

.3974 

.0094 

-.0109 

.0004 

-.0009 

4.4 

.3991 

.0074 

-.0086 

.0003 

-.0006 

4.6 

.4004 

.0059 

-.0068 

.0002 

-.0004 

4.8 

.4015 

.0047 

-.0054 

.0001 

-.0002 

5.0 

.4023 

.0037 

-.0042 

.0001 

-.0001 

5.4 

.4035 

.0024 

-.0026 

.0001 

-.0001 

5.8 

.4043 

.0015 

-.0016 

.0000 

.0000 

6.4 

.4049 

.0008 

-.0008 

.0000 

.0000 

7.0 

.4053 

.0005 

-.0004 

.0000 

.0000 

8.0 

.4056 

.0002 

-.0001 

.0000 

.0000 

9.0 

.4058 

.0001 

.0000 

.0000 

.0000 

10.0 

.4059 

.0001 

.0000 

.0000 

.0000 

11.0 

.4059 

.0001 

.0000 

.0000 

.0000 
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TABLE  I  -  FUNCTIONS  F  AND  H  AND  DERIVATIVES 
(f)  Prandtl  number,  10  (g)  Prandtl  number,  100 


n 

F 

F' 

F' * 

H 

•H* 

0.000 

0  .000 

0 .0000 

0.2517 

1.0000 

-2.191 

.025 

.0001 

.0060 

.2274 

.9452 

-2.191 

.050 

.0003 

.0114 

.2044 

.8905 

-2.188 

.075 

.0006 

.0162 

.1828 

.8359 

-2.180 

.100 

.0011 

.0205 

.1626 

.7815 

-2.166 

.125 

.0017 

.0244 

.1438 

.7276 

-2.144 

.150 

.0023 

.0277 

.1263 

.6744 

-2.113 

.175 

.0031 

.0307 

.1101 

.6221 

-2.071 

.200 

.0039 

.0332 

.0952 

.5709 

-2.018 

.225 

.0047 

.0355 

.0816 

.5213 

-1.954 

.250 

.0056 

.0373 

.0692 

.4733 

-1.880 

.275 

.0066 

.0389 

.0580 

.4273 

-1.786 

.300 

.0076 

.0402 

.0479 

.3836 

-1.704 

.325 

.0086 

.0413 

.0389 

.3422 

-1.604 

.350 

.0096 

.0422 

.0309 

.3034 

-1.498 

.375 

.0107 

.0429 

.0238 

.2674 

-1.388 

.400 

.0118 

.0434 

.0176 

.2341 

-1.276 

.425 

.0129 

.0438 

.0123 

.2036 

-1.163 

.450 

.0140 

.0440 

.0076 

.1759 

-1.052 

.475 

.0151 

.0442 

.0036 

.1510 

-.9430 

.500 

.0162 

.0442 

.0002 

.1287 

-.8393 

.525 

.0173 

.0442 

-.0026 

.1089 

-.7404 

.550 

.0184 

.0441 

-.0050 

.0916 

-.6478 

.575 

.0195 

.0439 

-.0070 

.0765 

-.5621 

.600 

.0206 

.0437 

-.0087 

.0634 

-.4837 

.65 

.0227 

.0432 

-.0111 

.0427 

-.3496 

.70 

.0249 

.0426 

-.0126 

.0280 

-.2440 

.75 

.02,70 

.0420 

-.0135 

.0178 

-.1657 

.80 

.0291 

.0413 

-.0140 

.0111 

-.1088 

.85 

.0311 

.0406 

-.0142 

.0067 

-.0693 

•  .90 

.0331 

.0399 

-.0142 

.0039 

-.0428 

.95 

.0351 

.0392 

-.0141 

.0022 

-.0256 

1.00 

.0371 

.0385 

-.0140 

.0012 

-.0149 

1.10 

.0408 

.0371 

-.0136 

.0004 

-.0046 

1.20 

.0445 

.0357 

-.0132 

.0001 

-.0013 

1.30 

.0480 

.0344 

-.0128 

.0000 

-.0003 

1.40 

.0514 

.0332 

-.0124 

.0000 

-.0001 

1.50 

.0546 

.0320 

-.0119 

.0000 

.0000 

1.60 

.0578 

.0308 

-.0116 

.0000 

.0000 

1.70 

.0608 

.0297 

-.0112 

.0000 

.0000 

1.80 

.0637 

.0286 

-.0108 

.0000 

.0000 

1.90 

.0665 

.0275 

-.0104 

.0000 

.0000 

2.0 

.0692 

.0265 

-.0101 

.0000 

.0000 

2.1 

.0718 

.0255 

-.0097 

.0000 

.0000 

2.2 

.0743 

.0245 

-.0094 

.0000 

.  .0000 

2.3 

.0767 

.0236 

-.0091 

.0000 

.0000 

2.4 

.0790 

.0227 

-.0088 

.0000 

.0000 

2.6 

.0834 

.0210 

-.0082 

.0000 

.0000 

2.8 

.0874 

.0195 

-.0076 

.0000 

.0000 

3.0 

.0912 

.0180 

-.0071 

.0000 

.0000 

n 

F 

F’ 

FT  r 

H 

H, 

0 

0 .0000 

0,0000 

0.4192 

1.0000 

-1.1694 

.1 

.0019 

.0371 

.3251 

.8831 

-1.1671 

.2 

.0071 

.0654 

.2428 

.7670 

-1.1521 

.3 

.0147 

.0861 

.1723 

.6534 

-1.1155 

.4 

.0241 

.1003 

.1134 

.5448 

-1.0526 

.5 

.0346 

.1091 

.0655 

.4437 

-.9640 

.6 

.0458 

.1137 

.0279 

.3527 

-.8545 

.7 

.0573 

.1150 

-.0011 

.2733 

-.7322 

.8 

.0687 

.1137 

-.0221 

.2064 

-.6061 

.9 

.0800 

.1107 

-.0367 

.1519 

-.4849 

1.0 

.0908 

.1066 

-.0462 

.1090 

-.3753 

1.1 

.1012 

.1016 

-.0518 

.0763 

-.2813 

1.2 

.1111 

.0963 

-.0545 

.0522 

-.2045 

1.3 

.1205 

.0908 

-.0552 

.0349 

-.1445 

1.4 

.1293 

.0853 

-.0544 

.0228 

-.0993 

1.5 

.1376, 

.0799 

-.0527 

.0146 

-.0665 

1.6 

.1453 

.0748 

-.0505 

.0092 

-.0435 

1.7 

.1525 

.0699 

-.0480 

.0056 

-.0278 

1.8 

.1593 

.0652 

-.0453 

.0034 

-.0174 

1.9 

.1656 

.0608 

-.0427 

.0020 

-.0107 

2.0 

.1714 

.0567 

-.0400 

.0012 

-.0065 

2.1 

.1769 

.0528 

-.0375 

.0007 

-.0038 

2.2 

.1820 

.0491 

-.0351 

.0004 

-.0022 

2.3 

.1868 

.0458 

-.0328 

.0002 

- . 0013 

2.4 

.1912 

.0426 

-.0306 

.0001 

-.0007 

2.5 

.1953 

.0396 

-.0286 

.0001 

-.0004 

2.6 

.1991 

.0369 

-.0267 

.0000 

-.0002 

2.7 

.2027 

.0343 

-.0249 

.0000 

-.0001 

2.8 

.2060 

.0319 

-.0232 

.0000 

-.0001 

2.9 

.2090 

.0297 

-.0216 

.0000 

.0000 

3.0 

.2119 

.0276 

-.0201 

.0000 

.0000 

3.1 

.2146 

.0256 

-.0187 

.0000 

.0000 

3.2 

.2170 

.0238 

-.0174 

.0000 

.0000 

3.3 

.2193 

.0221 

-.0162 

.0000 

.0000 

3.4 

.2215 

.0206 

-.0151 

.0000 

.0000 

3.6 

.2253 

.0178 

-.0131 

.0000 

.0000 

3.8 

.2286 

.0153 

-.0113 

.0000 

.0000 

4.0 

.2314 

.0132 

-.0098 

.0000 

.0000 

4.2 

.2339 

.0114 

-.0084 

.0000 

.0000 

4.4 

.2360 

.0098 

-.0073 

.0000 

.0000 

4.6 

.2379 

.0085 

-.0063 

.0000 

.0000 

4.8 

.2394 

.0073 

-.0054 

.0000 

.0000 

5.0 

.2408 

.0063 

-.0047 

.0000 

.0000 

5.4 

.2430 

.0047 

-.0035 

.0000 

.0000 

5.8 

.2446 

.0035 

-.0026 

.0000 

.0000 

6.2 

.2458 

.0026 

-.0019 

.0000 

.0000 

7.0 

.2474 

.0014 

-.0011 

.0000 

.0000 

8.0 

.2484 

.0007 

-.0005 

.0000 

.0000 

9.0 

.2489 

.0003 

-.0002 

.0000 

.0000 

10.0 

.2491 

.0002 

-.0001 

.0000 

.0000 
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FOE  VAEIOUS  PEAKDTL  BOMBERS  -  CONCLUDED 


(g)  Erandtl  number ,  100  -  Concluded 


(h)  Prandtl  number,  1000 


n 

F 

F’ 

F*  » 

H 

H* 

0. 

b.oooo 

0.0000 

0.1450 

1.0000 

-3.966 

.025 

.0000 

.0033 

.1212 

,9009 

-3.962 

.050 

.0002 

.0061 

.0999 

,8021 

-3.933 

.075 

.0003 

.0083 

.0811 

.7046 

-3.861 

.100 

.0006 

.0102 

.0647 

.6096 

-3.731 

.  .125 

.0008 

.0116 

.0506 

.5186 

-3.538 

.150 

.0012 

.0127 

.0387 

.4332 

-3 .283 

.175 

.0015 

.0135 

.0289 

.3549 

-2.975 

.200 

.0018 

.0142 

.0209 

.2847 

-2.628 

.225 

.0022 

.0146 

.0146 

.2236 

-2.261 

.250 

.0026 

1,0149 

.0096 

.1717 

-1.893 

.275 

.0029 

.0151 

.0059 

.1288 

-1.541 

.300 

.0033 

.0152 

,0032 

.0944 

-1.220 

.325 

.0037 

.0153 

.0012 

.0675 

-.9381 

.350 

.0041 

.0153 

-.0003 

.0471 

-.7012 

.375 

.0045 

.0152 

-.0012 

.0321 

-.5093 

.400 

.0048 

.0152 

-.0019 

.0213 

-.3596 

.425 

.0052 

.0152 

-.0023 

.0138 

-.2467 

.450 

.0056 

.0151 

-.0026 

.0087 

-.1645 

.475 

.0060 

.0150 

-.0027 

.0053 

-.1066 

.500 

.0063 

.0150 

-.0028 

.0032 

-.0672 

.525 

.0067 

.0149 

-.0029 

.0019 

-.0412 

.550 

.0071 

.0148 

-.0029 

.0011 

-.0245 

.575 

.0075 

.0147 

-.0029 

.0006 

-.0142 

.600 

.0078 

.0147 

-.0029 

.0003 

-.0080 

.625 

.0082 

.0146 

-.0029 

.0002 

-.0044 

.800 

.0107 

.0141- 

-.0028 

.0000 

.0000 

1.000 

.0135 

.0136 

-.0027 

.0000 

.0000 

1.40 

.0187 

.0125 

-.0025 

.0000 

.0000 

1.80 

.0235 

.0115 

-.0023 

.0000 

.0000 

2.20 

.0279 

.0106 

-.0022 

.0000 

.0000 

2.60 

.0320 

.0098 

-.0020 

.0000 

.0000 

3.0 

.0358 

.0090 

-.0019 

.0000 

.0000 

3.6 

.0409 

.0080 

-.0017 

.0000 

.0000 

4.2 

.0454 

.0070 

-.0015 

.0000 

.0000 

5.0 

.0505 

.0060 

-.0012 

.0000 

,0000 

5.8 

.0549 

.0050 

-.0011 

.0000 

.0000 

7.0 

.0603 

.0039 

-.0008 

.0000 

.0000 

8.0 

.0638 

.0032 

-.0007 

.0000 

.0000 

10.0 

.0691 

.0022 

-.0004 

.0000 

.0000 

12.0 

.0727 

.0015 

-.0003 

.0000 

.0000 

14.0 

.0752 

.0011 

-.0002 

.0000 

.0000 

16.0 

.0771 

.0008 

-.0001 

.0000 

.0000 

18.0 

.0786 

.0007 

.0000 

.0000 

.0000 

20.0 

.0798 

.0006 

.0000 

.0000 

.0000 

22.0 

.0809 

.0005 

.0000 

.0000 

.0000 

23.6 

.0816 

.0005 

.0000 

.0000 

.0000 

(U)H 


Figure  2.  -  Dimensionless  temperature  distributions  for 
various  Prandtl  numbers. 
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Figure  3.  -  Comparison  of  small  plate  experimental  and  theoretical 
velocity  distributions  for  Prandtl  number  of  0.72. 
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Figure  5.  -  Comparison  of  large  plate  experimental  and  theoretical 
velocity  distributions  for  Prandtl  number  of  0.72. 
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temperature  distributions  for  Prandtl  number  of-  0.72. 
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NACA-Langley  -  2-21-52  -  1000 
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